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APPLICATION OF RADIAL-EQDIT T K RIDM CONDITION TO AXIAL-FLOW 
COMPRESSOR AND TOEBINE DESIGN 
By Chung-Eua Wu and Lincoln Wolfenstein 


SUMMARY 

Basic general equations governing the three-dimensional compres- 
sible flow of gas through a compressor. or a turbine are given in 
terms of velocity components, total enthalpy, and entropy. These 
equations are used to determine the radial motion of gas through an 
axial-flow compressor or a turbine and the corresponding effect on 
the radial variations of the state of gas between successive blade 
rows in the case of steady, axially symmetrical flow. The aspect 
ratio of the blade row is found to be an important factor in the cal- 
culation when the effect of radial motion is included. The usual 
method, which negleots the effect of radial motion, is shown to be 
good only for the limiting case of zero blade-rcw aspect ratio, that 
is for the case where the axial length of the blade row is much 
larger than the radial length of the blade row. A sinusoidal radial- 
flow path is found to give the effeot of radial motion on the radial 
variation of gas state between blade rows as small as likely without 
any discontinuity in the curvature of the streamline and is sug- 
gested for use in design calculations. 

The equations are applied to investigate the maximum compatible 
number of the radial variations of the gas properties between 
successive blade rows that a designer is free to specify. The vari- 
ous ways of taking up these degrees of freedom and the different 
types of design obtained are discussed. A general procedure is given 
to calculate the characteristics of a compressor or turbine of any 
given type of design, taking into acocunt the effect of radial 
motion of gas. Numerical calculations made for two types of com- 
pressor and one type of turbine show that even in the case of non- 
tapered passage, there is appreciable radial motion and that the 
corresponding effeots are of significant magnitude ana , should be 
taken into acoounb in design. 
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INTROKJCTION 


The design of a compressor or a turbine hereinafter called a 
turhccoachlne may he divided into two phases. The first phase con- 
cerns the type of design to he used or the determination of the 
most desirable variations of velocity and therm odynamio properties 
of the gas in planes normal to the axis of the machine between 
successive blade rows. The second phase concerns the design of 
blades that will give the desired variations of velocity and, other 
properties of gas in these planes. In the first phase, the condi- 
tion of radial equilibrium, that is, the radial component of the 
equation of motion, must be used. The flow of gas in a turbo- 
machine is curvilinear; it is curved not only by the whirling 
motion of gas, but also by the radial motion of the gas. The equa- 
tion of motion then specifies the radial pressure gradient required 
to provide the centripetal foroe to maintain the curved flow. 


In figure 1(a), a curved stream surface over four stages of a 
multistage turbamachine is shown and in figures 1(b) and 1(c) are 
shown its intersections with planes normal to and containing the 
axis of the machine, respectively. The radial pressure gradient 
d ue to the whirling motion of gas is always positive, whereas that 
due to the radial motion of gas may be either positive or negative 
depending on whether the curvature caused by the motion is either 
Inward or outward from the axis of the machine at the point of 
consideration. If the gas lying on a normal plane between two 
blade rows, such as station 1 in figure 2(a), is considered, the 
equation of motion gives (neglecting small terms) 


p Ir 



dV 

r 

"ST 


- - v. 


2 

2 d r 

: ^ 2 
dz 


where p arid p are the pressure and the de ns ity of the gas 
particle, respectively, r and z are the radial distance from 
and the distance along the axis of the gas particle, respectively, 
and v r , Vq, and V z are the radial, tangential, and axial 
components of the velocity of the gas particle, respectively. 

From this equation it iB seen that the effect of the two motions 
on the radial pressure gradient is quite similar. The effects of 
the two motions are both proportional to the product of the square 
of the velocity and the curvature (l/r and d 2 r/dz , respec- 
tively) involved. Even when the radial motion involved is small, 
if the axial velocity is high and the blade-row aspect ratio is 
large, the second term in the equation is of comparable magnitude 
to the first and should be included in the calculation. In the 
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past, however, the radial motion and its effect were usually 
neglected (references 1 to 3). (Such methods of calculation, will 
hereinafter he referred to as the "sinplif led- radial- equilibrium 
approximation.") In reference 1, a general form of the Euler 
equation is presented hut is not used later on in the through-flow 
analysis because the simplified- radial-equilibrium consideration 
may he expected to hold sufficiently far downstream of a single row 
of blades where the radial velocity of the gas becomes small. Eor 
this reason, this approximation is adequate to predict most of the 
experimental results available on the radial variation of the gas 
state far downstream of a single row of blades (references 1 
and 3) . No satisfactory data or theory exist for the. general case 
of a blade row within a closely spaced series of blade rows, as in 
a multistage compressor or turbine where the effect of radial 
motion may be quite significant. In the present paper, the effect 
of the radial motion of gas is considered in applying the radial- 
equilibrium condition to turbomachine design. 

In the ANALYSIS, the general equations governing the flow of 
gas in axial-flow turbcmachines are developed primarily for the 
case of steady axially symmetrical flow, which corresponds to the 
limiting case of an infinite number of blades. The blades are 
replaced in the calculation by an appropriate force field. A 
method of solution involving the use of the basic equations in 
finite-difference form is discussed. Expressions are developed 
for the use of a large number of successive axial stations and for 
the use of three stations for a stage in which an appropriate 
radial-flow path is assumed. 

The equations are applied to investigate the maximum number of 
radial variations of the velocities and other properties of the gas 
that a designer is free to specify in any particular design. The 
various ways of using these degrees of freedom and the different 
types of design obtained are discussed. 

In the last part of this report, a general procedure of cal- 
culation for any type of design is given. It Is used to calculate 
two types of compressor and one type of turbine with the purpose 
of investigating the effect of radial motion on design calculations 

This investigation, conducted at the HACA Cleve land laboratory 
was completed in April 1948. 
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SYMBOLS 

The following symbols are used in this report: 

aspect ratio of blade row, 

velocity of sound 
lift coefficient 
specific heat of gas at constant pressure 
specific heat of gas at constant volume 

differentiation with respect to time following motion of gas 

force acting on gas particles by blades per unit mass of gas 

radial component of F 

tangential component of F 

axial component of F 

form of radial-flow path 

mass flow per unit flow area perpendicular to axis of 
turbomachine 

form of radial-displacement distribution 
total enthalpy per unit mass of gas, 
enthalpy per unit mass of gas, ^u + 
constant 

axial length of blade row (fig. 2) 

Mach number of gas 
mass of gas 
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n 


P 

Q 


R 

r 


r ' = 


polytropic exponent of actual expansion or compression 
process of gas 

static pressure 

external heat transferred to gas particle along its path of 
motion per unit mass per unit time 

gas constant 

radius measured from axis of turboma chine 

mean radius of radial-f lov path (fig. 2(a)) 
r 


s entropy per unit mass of gas 

s.r.e. simplified- radial-equilibrium approximation 
T absolute stream temperature of gas 

t time 

U magnitude of U 

U vector velocity of rotor blades at radius r 

u internal energy per unit mass of gas with 0° absolute as 

base temperature 

V magnitude of V 

Y absolute vector velocity of gas 

Y r radial component of Y 

Y_ tangential component of Y 

0 

Y„ axial component of ? 

<£t 

W magnitude of W 

W vector velocity of gas relative to rotor blade 
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tangential component of W 

maximum radial displacement over "blade height 
distance along axis of turbama chine 

angle "between relative velocity of gas entering rotor and 
axis of turboma chine 

angle "between absolute velocity of gas entering stator and 
axis of turbomachine 

ratio of specific heats, (cp/cy) 

radial displacement across rotor, (r 2 - r-^) 

^ 2 - 5 1 \ 

dimensionless turning, I — = I 

\ r l u l,t/ 


small-stage or polytropic efficiency 

angular distance measured from some fixed radial line 

viscosity of gas 

kinematic viscosity of gas, (jx/p) 

angular momentum about z axis per unit mass of gas, 


(rVj 


p mass density 

O blade solidity 

<ji dissipation of energy due to viscosity per unit volume of 

gas per unit time 

cp function 

00 angular velocity of rotor 

Subscripts: 

1 in front of rotor 


2 


behind rotor and in front of stator 
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3 behind stator and in front of next rotor 

4 behind next rotor 

0 satisfying continuity equation 

e satisfying radial- equilibrium and total- enthalpy equations 

h at hub 

1 any station 

J station short distance downstream of station i 

k any station in front of rotor 

L limiting value 

7 any station in front of stator 

m used with r to Indicate mean radius 

n used with r to indicate radius where maximum radial 

displacement occurs 

s simplified-radial-equilibrium approximation 

t at tip 


ANALYSIS . 

Basic Equations 

The state of gas in three-dimensional motion is completely 
specified by its absolute vector velocity Y, or its three com- 
ponents Y r , V q, and Y z referred to cylindrical coordinates r, 
6, and z, respectively, and two thermodynamic properties. These 
last two properties are usually chosen as the static pressure p 
and the static density p, but herein it is more convenient to use 
the total enthalpy per unit mass H and entropy per unit mass s, 
which are defined by 


H 



( 1 ) 
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and 


T ds = du + p dp" 1 (2) 

For the range of temperature and pressure encountered in compressors 
and turbines/ p, p, and T are accurately related by the fol- 
lowing equation of state 

p = pET (3) 

The Navi er- Stokes equation of motion for a real fluid is 
given, in vector form, by 

B? = F- ivp+l± [v 2 V + 1 V(V*Y)1 (4) 

Dt PPL 3 J 

The energy equation for a real fluid is given by (refer- 
ences 4 and 5) 

£* + p 5P^: = Q + t (5) 

Dt Dt P 

The continuity equation can be written as 

(py) = 0 ( 6 ) 

or 

V* V + ^ log Q p = 0 (6a) 


When the preceding equations are combined, the following four 
general equations are obtained (for derivation, see appe n d ix A): 


VH = F + TVs + V X (V X V) 




™ = q + ^ + 1 |£ 
Dt p p dt 


+ V* < 


F + 


v Jv 2 ? + | 


v(v-v) 


(?) 

( 8 ) 


Ds 

Dt 



+ E $- 
P 


o 
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7 - ?+ ^IS 1 °Se T -i(l )= 0 < 10 > 

Equation (7) gives the gradient of total enthalpy in terms of "blade 
force, viscous forces, velocity, and other properties of the gas. 
This vector equation gives three scalar equations in three dimen- 
sions. Equation (8) gives the rate of change of total enthalpy of 
gas along a streamline in terms of rate of beat addition, rate of 
work done by blade and viscous forces, and so forth. Equation (9) 
gives the rate of change of entropy along a streamline in terms of 
rate of external heat addition and of dissipation of energy due to 
viscosity. Equation (10) gives the continuity relation in terms 
of velocity, temperature, and entropy of gas. It should be noted 
that of the six scalar equations given by equations (7) to (10) 
only five are Independent, because they are derived from equa- 
tions (4) to (6), among which equation (4) yields three scalar 
equations, giving a total of five Independent equations. (The 
energy relation. is used in the derivation of both equations (8) 
and (9).) 


Steady Axially Symmetric Flow- 

Axial symmetry, that is symmetry about the axis of rotation, 
can be assumed to exist sufficiently far downstream of any blade 
row and is true everywhere for the limiting case of an infinite 
number of blades. This assumption is usually made to simplify 
the analysis. The blades are then replaced by a volume distribu- 
tion of forces, the magnitude of which iB obtained by maintaining 
constant the product of the number of blades and the resultant 
force at any point on the blade as the number of blades is 
increased. (The resultant force is the difference between the 
forces on the two sides of the blade.) That is, the resultant 
force acting on the gas by a blade element at any radius is con- 
sidered to be evenly distributed over the stream sheet between 
two blades at that radius. Reference 1 shows that for incompress- 
ible and frictionless flow the value thus obtained gives an aver- 
age value with respect to the coordinate 6 . Because the number 
of blades is usually large, this assumption is considered reason- 
able for steady operation of turbomachines, and particularly for 
this investigation. Thus all partial derivatives of gas prop- 
erties with respect to time and angular position 6 will be 
taken as equal to zero. The state of gas is considered as a 
function of r and z only, that is, the problem is reduced to 
a two-dimensional treatment. With this simplification and by 
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transfoiming V 2 Y into £v(V*y) - 7x (V X Y)], there are 

obtained the radial, tangential, and axial component a of equa- 
tion (7): 


Sh _ -p . m Ss Y e s ( rT e) 

to ~ r &? “ ~*r 


Sy z 5v r 
+ Yz !& “ Vz aT 


+ 7 


4 d 2 (rV r ) 4 5(rY r ) ± b\ 

‘ 3r 2 “ 5F“ + 3 5rSz + Sz 2 


3r Sr 2 


(7a) 


T r d(rY e ) bV Q 

0 cF e - r Sr “ Yz Bz + U 


1 d 2 (rV e ) ! 5(rY e ) d 2 ^ 


r Sr 2 r- 




and 


(7b) 


SH _ _ ^ 5s _ gr _ jgz 

" F z + T 5z + Vr 5z “ Y * 5r + T 0 ^ 


5z 


+ v 


5 2 (rY r ) 4 5^ bHj. d^Y 


[JV 4 

3r 5r5z 3 ^ 2 


drdz ^r 2 


(7c) 


In the radial equilibrium equation (7a), the relative impor- 
tance of the various terms depends mainly on the type of design. 

The radial force exerted by the blade depends an the twist and the 
taper of Idle blade in the radial direction. The radial variation 
of entropy depends on the radial variations of external heat sub- 
traction and of dissipation of energy due to' viscous effect. In 
ordinary compressors and turbines, the amount of heat transfer is 
small per unit mass of gas flow, and its effect on the radial var- 
iation of entropy is negligible. But in turbines with blade cool- 
ing, the heat transfer may be quite large and it may give a signif- 
icant radial variation of entropy, which must be taken into consid- 
eration. The increase of entropy due to dissipation of energy by 
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bJ 

a 

O) 


viscosity is small and may "be assumed to be the same at all radii 
over the main portion of gas outside the boundary layers at the 
rotor drum and the inner wall of the outer casing. The third term 
on the right-hand side of the equation is equal to zero for free- 
vortex type of design, in whioh rVg is set constant at all radii, 
gnfl has a nonzero value in other types of design. Similarly, the 
fourth term may have a nonzero value if constant axial velocity at 
all radii is not specified in design. 

The fifth term represents the effect of radial motion on the 
radial-equilibrium condition. In the gap between two blade rows 
it is proportional to the product of the square of axial velocity 
and the curvature of the flow path in the axial plane and is very 
similar to the effect of whirling motion represented by V 0 2 /r. 

If the curvature caused by radial motion is positive, the effect 
is to decrease the radial gradient of total enthalpy or pressure 
caused by the whirling motion of the gas. If the curvature is nega- 
tive, the effect is to intensify this radial gradient. For design 
with large blade height, short axial blade length, high whirl, 
e/nfl high axial velocities, it will be shown that the effect of 
radial motion is large and should be considered in the calculations. 


The last four terms in the bracket are usually of the same 
order of magnitude as the preceding velocity derivatives and 
because they are multiplied by u, which is much smaller than 
other multipliers (V z or Vg), the whole product is much smaller 
than the other terms in the equation and may be neglected. The 
same argument applies to similar terms in equations (7b) and (7c). 
Hence equations (7a), (7b), and (7c) may be simplified to 


a n d 


SH 

Sr 


F^ + 


T 

Sr r Sr z Sr 


- 


dv x 

Sz 


0 



T r S(rVg) SVg 
r Sr ^ z Sz 


SH 

Sz 


SB = f„ + T 


Ss 

Sz 


+ V r 


Sv 

r 

Sz 


- V„ 


Sr 


+ V, 


Sz 


(7d) 

(7e) 


(7f) 


* 
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Among the terms on the right-hand side of equation (8) , other 
than the first term (the importance of which depends on whether 
there is blade cooling), Y*F is the predominating term in passing 
through the rotor. Whereas in passing through the stator, all the 
terms on the right-hand side are of the same order of magnitude. 
However, by the use of the equation of motion, the energy equation 
for steady flow, and the assumption that the heat generated from 
the frictional work remains in the stream, DH/Dt can be expressed 
for steady axially symmetrical flow in a more useful form in terms 


of only Q, tt) , 


and 


D(rV e ) 

Dt 


(see appendix A): 


EH 

Dt 


Q +co 


D(rV e ) 

Dt 


(8a) 


In equation (9), the two terms are of comparable magnitudes 
in ordinary compressors or turbines. With blade cooling, the first 
term is predominate. For the axially symmetrical flow, it reduces 



For axially symmetrical flow, equation (10) reduces to the 
following 


1 d(rY r ) dY z p 
r Sr + & + 7-1 


^ loge T + Y z 




(10a) 


Among the preceding six equations, equations (7d) to (10a), 
only five are independent, because both equations (8a) and (9a) 
represent the energy relation. For nonviscous flow, there is an 
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additional relation in that the blade force ia normal to the sur- 
face of the blade and consequently is perpendicular to the relative 
velocity of the gas, that is, 

F • (V - U) = 0 (11) 


Then equation (8a) can be obtained by using the motion and energy 


equations and equation (11). 
either equation (8a) or (ll) 
additional relation, giving, 
(equations (7d), (7e), (7f), 


For such flow, equation (9a) reduces to simply 


(See appendix A. ) Consequently, 
can be considered as representing this 
in all, six independent equations 
(9a), (10a), and either (8a) or (ll)). 


Ds = Q, 
Dt T 


(9b) 


For viscous flow, equation (ll) is no longer true because the 
force exerted by the blade is now inclined from the direction per- 
pendicular to the relative velocity of the gas by an amount to be 
determined from the shearing stresses of the gas adjacent to the 
blade surface. Under the present consideration of axially 
symmetrical flow it is found desirable to retain equation (8a) as 
the energy equation and to obtain the entropy change from H, V, 
and the polytropic exponent n of the actual compression or 
expansion process: 


p p -n = constant (12) 

Then the rate of change of entropy along the streamline is (see 
appendix A) 

ll = E Tn-ifir-iy ( v * I? loge T + v z 'k loge *) (13) 

In equation (13 ) n is considered as known. In a given machine, 
n may be obtained directly from measured pressure and temperature 
data. In a new design, n may be obtained from the assumed poly- 
tropic efficiency used in design calculations: 

For the compressor 




2=1 


_Z_ 

n- 


or 


n = — 
1 


1 

. 12= 1 


n 7 


n 
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For the turbine 


n-1 



7 


n 


E= 


1 • 



Because the change in s is usually small compared with the ch an ges 
in H and Y, the preceding method of determining s seems to he 
adequate to account for viscous effect in calculating the pressure 
and, density change along the streamline for the present problem. 


Method of Solution 

If equations (l), (2), and (3) are considered as the relations 
to express p, p, and T in terms of H and s, then equa- 
tions (7d), (7e-), (7f), (8a), (9b), and (10a) are six independent 
equations in the case of nonviscous flow, and equations (7d), (7e), 
(7f), (8a), (10a), and (13) are six independent ^equations in the 
case of viscous flow both involving eight variables v r> T z> 

H, s, F r , F e , and F z . In the direct problem, in which the 
shape of the blade profile, the shape of the inner and outer wall 
of the gas passage, the rotor speed, the power input, and the 
entering and exit conditions of the gas are given. It is theoret- 
ically possible to determine the variations of these quantities 
throughout the machine. In the inverse problem, in which the 
desirable variations of two of the gas properties are pre- 
scribed, It Is also theoretically possible to determine the varia- 
tions of the other properties of the gas and the blade force 
necessary to achieve the prescribed variation of gas conditions. 
However, it seems that no general analytical solution of these 
equations is possible in either problem. Two numerical methods of 
solution are therefore suggested. In the first method, the 
previous equations are expressed in finite-difference form and 
applied to successive axial stations that are a short distance 
apart. In the second method, only three stations are used for 
each stage in which an appropriate radial-flow path is assumed. 

Method of finite difference for successive axial stations. 

At each station, if tYq is denoted by £, equations (7d), (7e), 
and (7f ) can be written as 
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av r \ 
" -SV 



(7g) 

(7h) 


and 



Between any two successive stations i and J that are a 
short distance apart (fig. 3), the change in total enthalpy is • 
given "by equation (8a) 


Hj(rj) - H^ri) «» (o |jj(rj) - + P J Q dt 

vJ-ti 


(8h) 


where (r) indicates that the gas properties at a particular sta- 
tion are a function of the radial position of the gas particle in 
that station; (It should he noted that due to radial motion the 
radial position of a gas particle at any station J is different 
from its radial position at the previous station i.) In passing 
through the stator, co = 0 and there is no change of H along 
each streamline except from heat effect. 

The entropy change between the two stations is obtained from 
equation (13) (see appendix A); 


s j( r j) - - s log e 



(13a) 


Instead of integrating equation (10a), the continuity relation 
between the two stations is readily obtained by equating the mass 
flow at the two stations. 


Gj rj drj = Gi r* dr* 


(14) 
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By expressing G in H, Y, and s, 
appendix A) : 


equation (14) "becomes (see 





. S r., ax } = V 2)I 



dr. 


(14a) 


Equations (7g), (7h), (7i), (8b), (13a), and (14a) are now six 
independent equations relating tbe gas properties and blade forces 
at tbe two stations. In these equations, the heat-transfer term is 
negligible in ordinary turbomachines and can be estimated in the 
case of cooled blades; the temperature I is a known function of H 
and Y; and r^ is obtained from r^ and Y rj Hence there 
are only eight unknowns in Hj, sj, Y r# j, Y^j, Y Z>J , F r ^ j, 

F q j, and F^ j (those at the first station i are considered 

as known) . For a given blade operating at a given speed, two 
additional relations are known among idle velocity components from 
the tangent plane to the blade surface at the point (rj, z j) 
because the flow of gas has to conform to the shape of the blades. 
For a new design, the desirable variation of two of the eight 
variables (usually one of the two is Y g ) can be specified and the 
remaining five determined from the preceding equations. It Is to 
be noted, however, that, in practice, the radial blade force 1 b 
not essentially an independent variable to be specified by the 
designer, but Is mainly determined by the actual construction of 
the blade to meet the aerodynamic as well as the mechanical- 
strength requirements at different radii. The designer has to see 
to It, of course, that the radial force to be obtained from the 
blade actually constructed Is consistent with that used in or 
obtained from the design calculations. The procedure of calcula- 
tion for this step-by-step method varies with the type of design, 
the conditions given, or the two gas properties prescribed. The 
calculation Is quite laborious and seems to be Justified only in 
the process of actual design. 


In order to obtain an over-all picture of the radial motion 
in a turbomachine and its effect on design considerations, the fol- 
lowing method of using only three stations for each stage with pre- 
scribed radial-flow path may be used. 
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Method of prescribed radial-flow path. - In a turbomachine, 
the radial motion of the gas is caused by three factors: 

(1) Tapering of the annular passage either at the inner or 
outer -Ball gives the flow a radial displacement across the stage, 
which is, of course, greatest in the immediate neighborhood of the 
tapered surface. 

(2) Even with a nontapered passage, a radial displacement 
across the stage may be called for because of a variation in the 
distribution of specific mass flow over the blade height. 

(3) Even if there is no radial displacement across the stage 
(that is, the same particle occupies the same radial position at 
the first station of each successive stage), there will, in 
general, be radial displacement of flow within the stage. This 
radial flow will then be oscillatory in nature, a radial displace- 
ment in the rotor being followed by an equal and opposite radial 
displacement in the stator. This radial flow arises because of 
the difference between the radial variation of the specific mass 
flow within the stage and those at the entrance and exit stations 
of the stage. (This radial displacement can only be avoided by 
specifying zero or the same radial variation of specific mass flow 
at all stations of the stage in the design.) 

In general, the radial flow of gas therefore consists of a 
gradual, generally monotone, radial motion due to factors (l) and 
(2) with an oscillatory motion of period equal to the stage length 
due to factor (3) superimposed on it. The radial flow caused by 
these three factors will be similar to that shown in figure 1. 

The effect of the radial motion an the calculations arises chiefly 
through the term dV r /<3z in the radial-equilibrium equation (7g). 
This term is expected to be significant mainly because of the 
oscillatory motion, which may require significant changes in Y r 
within a single row of blades. The case of oscillatory motion 
within a stage with no over-all radial displacement across the 
stage will therefore be considered first. That is, the gas passage 
is nontapered the radial distribution of gas properties at the 
entrance and exit stations of the stage is the same. 

Because there is no blade force acting on the gas and there is 
little time available for the gas to mix, the gas flowing through 
the gap between two blades is under a constant pressure gradient 
and consequently tends to move with the same curvature it acquires 
while leaving the first blade. For straight passages, the maximum 
and minimum points of the radial-flow path are likely to be 
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somewhere near the middle of the gap. (The intersecting curve of a 
stream surface with an axial plane is herein referred to as "the 
radial-flow path." Because of axial symmetry, the radial-flow 
path is the same in any axial plane.) The stations between blade 
rows are most conveniently chosen at these points. The stations 
in front of the rotor, between the rotor and the stator, and behind 
the stator are denoted by subscripts 1, 2, and 3, respectively. 

(See fig. 2(a).) If r Q and L represent the mean radial distance 
of the flow path and, the axial length of the blade row, respec- 
tively, then the radial distance of the gas particle at position z 
is given by 



at stations 1, 2, and, 3 

! - o, i, 2 

f (0) = f (2) = 1, f(l) = -1 1 (3,6) 

f ’(0) = f *(1) = f '( 2) =0 J 

where f is a function giving the form of the radial-flow path and 
the prime indicates differentiation with respect to z/L, It fol- 
lows that 


-rr dr 



f ' 



at station 1, z = 0, 


and 






(17) 


(18) 


inasmuch as bY z /bz Is practically zero in passing through the gap. 
Similarly, at station 2, z = L, 
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and 



0 




IT 


(18a) 


Because f" (z/L) determines bV r /bz or the effect of radial 
motion on the ra dial-eguilrbrium condition, it is desirable that 
it vary continuously; this condition together with those of equa- 
tion (16) suggests 


f 



«= COS It 



Then 


f"(0) = - it 2 , f"(l) = n 2 

and equations (15), (18), and (18a) become, respectively, 

r - r Q = - 2 g 1 cos it | (19) 


*V) r 2 ■* r l It 2 T 

i 2 l 2 2,1 


( 20 ) 



r 2 - r l it 2 v 

2 ,2 z > 2 


(20a) 


For the sinusoidal form of f(z/L), the maximum absolute 
value of f"(z/L) occurs at z = 0, L, and 2L, and is equal 
to it 2 . Even if f "(z/L) is assumed constant between z = 0 and 
z = l/ 2, thus minimizing the maximum absolute value of f” in 
the interval, the absolute value of f" equals 8. This assumption, 
however, necessitates a discontinuity in f" at z = L/2. The 
values of it 2 for the absolute values of f"(0) and f"(l) can be 
therefore considered as small as is likely. The smooth variation 
of f"(z/L) and the minimization of the absolute value of f "(z/L) 
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at the stations are reasonable assumptions provided that the radial 
force exerted by the blade remains relatively small. In such case, 
the sinusoidal curve is believed to represent the major harmonic of 
the actual radial-flow path, and the major effect of the radial 
motion may be obtained through the use of this simple curve. 


The radial- equilibrium equation (7g) may be written in terms 
of (r 2 - r-j_) by use of equation (20) 


m i m as i 1 . _ av z,i , lU l , » 2 

dr-L ~ Ti dr ± + ^2 ^ dr* + Y z,i dr^ + 2 ^ r 2 " r l^ ^2 T z,i 

( 21 ) 


ds., 






where 


1, 2, 3 


When r^ is replaced by the dimensionless variable r^' 


1 * 4 - - 




the dependence of this equation on the blade-row aspect ratio A is 
seen to be 


dHj 

drj 


r = Ti 


dsj 


t. 


k 


te ±' ri , 2 r t - % ir i' V r t - % 


+ t s ,i + C-U 1 | ( V - *1’) 

This form of the radial-equilibrium equation is seen to contain a 
term directly proportional to the radial displacement, to the 
square of the axial velocity, and to the square of the blade-row 
aspect ratio. If the blade-row aspect ratio is large or the axial 
velocity is high, the effect of radial motion may be large even 
though there is only a small amount of radial displacement across 
the blade row. 

This method is readily extended to the case where there is an 
over-all radial displacement across the stage due to tapering of 
the passage or due to variation in the design from stage to stage. 
In figure 2(b), the radial position of a gas particle originally 
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at r^ in station 1 is at r^ in station 3. For the oscillatory 
motion required within the stage, n °t generally equal to 

■| (r^ + r 3 ) . For the same reason stated in the previous case, it 

is desirable to have the radial-flow path consisting of a sinusoidal 
curve superimposed on the line passing through (z-^, r^) and 
(z 3, r 3 ); that is 


r , , r 5 - r l z , 1 f r l + y 5 w- * Z 

r-r l + 2 L + 2 V r 2 ” ° 1 “ 008 


( 22 ) 


Then 


y = — y s y 
v r cLz z z 


r 3 ' *1 if r l + r 3\ * 


2L + 2 l r 2 " 2 ) L Sin L 


JCZ 


(23) 


and 


-ffl 


M di( 

(*■1) 2 1 r 2 ” o / o F, 


2 v z,i 


inasmuch as BV z /dz is practically zero in passing through the gap. 
With this value of <5V v /bz, the radial- equilibrium equation (7g) 
becomes 



4si t i 

T-» . + p j_ + V . 

^i r ± 2 ^i z ^ 


dY z,l , _ 1 / fl^j3\ * 2 

dr ± 2 2 ) 2 


'z,i 


(24) 

where 

i = 1, 2, 3 

This equation is similar to equation (21). (If r^ = r-j_, it 
reduces to equation (21).) A similar equation in dimensionless r i * 
can also be obtained for this case by dividing r^ by (r^ - 1^)3.. 

At any point (z, r) within the blade region, the magnitude 
of F r consistent with this sinusoidal radial-flow path is obtained 
from equations (7g) and (23): 
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SH m ds i ^ 

35“ T 3?“33F" T z“5F + 




"2L“ z 


r 3 - r l 


r, + r- 


+ f l r 2 - 


f T i (> i °° s f + sto 


JtZ ^z 
L "ST 


(25) 


Degrees of Freedom in Design 

Before starting to design a turbomachine it is necessary to 
decide how the velocity and other properties of the gas should vary 
radially at successive axial stations to give the best type of 
design for the particular application. It is therefore necessary 
for the designer to know before hand the maximum compatible number 
of such variations that can be specified. 

In the discussion entitled "Method of finite difference for 
successive axial stations" it is seen that theoretically radial 
variation of two gas properties can be specified by the designer 
at each axial station within the blade region with the blade forces 
determined accordingly. In the free space between blade rows, the 
properties of gas remain constant along the streamline and no 
arbitrary change can be specified. In current design practice, 
the usual procedure is to specify the desirable gas conditions only 
at stations between blade rows (as well as before the first and 
after the last blade row) and then either to select some standard 
blade sections or to design them on the basis of prescribed veloc- 
ity distribution to achieve the desired change of gas state across 
the blade row. In either way, one degree of freedom at the second 
such station is taken up in going from the first station to the 
second station. The designer is therefore free to specify only one 
condition at each station throughout, the machine, with the exception 
of one station, usually the entrance station to the first stage. 

'Ways of specifying degrees of freedom. - In this discussion 
the following two f undam ental equations will be used: 

In the stations between blade rows, if the radial variation of 
entropy is negligible, the radial-equilibrium equation (7g) reduces 
to 

£% ii_£k v V (***) 

dr ± " 2 dr ± + v z,i dr-L “ v z,i \TzJ 1 


( 26 ) 
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Equation (8b) is applied to three successive stations of a stage 
for the case where' there is negligible heat transfer and is differ- 
entiated with respect to r-j_ 


ffi ffa SV fs ffs 

di^ + W \^dr 2 dr 1 " dr 1/ / “ dr 2 dr^ ” dr 3 dt^ 


(27) 


The following ways of taking up the degrees of freedom at 
these stations between successive blade rows sire discussed. 


(1) Constant work per unit mass of gas flow over the blade 
height. This condition is usually specified in the design of a 
turbomachine . It relates £ after the rotor to its value before 
the rotor in the following manner 

^2 (*2-) = h < r l> + ^t 5 t u l,t ( 2Q ) 


^ 2^2 ^1 

to 2 to l = to l 


(28a) 


where r-^ t 5^. U-^ . is equal to (£ g “ £q) & t tlie blade tip and 
is also equal to , (| g - | at other radii. 


Constant work over the blade height gives constant total- 
enthalpy change over the blade height. If the velocity at the 
exit of a stage is equal to that at the entrance, this condition 
also gives constant static- enthalpy cha n ge over the blade height. 


Under the condition of constant work, equation (27) reduces 
to 


^ = ^2 = f*3^ 

dr, ■ dr„ dr, dr„ dr, 

1 2 1 3 1 


(27a) 


(2) Constant total enthalpy over the blade height: 


SH 

3? = 


0 


(29) 
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This condition holds for the first stage of a compressor and will 
hold for all succeeding stages if constant work per unit mass over 
the blade height is employed. If a nonzero value of dH]_/dr-|_ is 
desired, an initial preparatory stage must be specially designed to 
obtain this value. In the last stage, however, it is usually 
desirable that SH/c)r be nearly zero. 

(3) Free-vortex-type distribution of tangential velocity: 

§ - 0 (30) 

or 


e i - r i T e,i = % 


(30a) 


This condition is commonly used in turbines and compressors. For 
incompressible flow, in addition to this condition, constant total 
enthalpy and constant axial velocity over the blade height can be 
specified. But for compressible flow, radial motion exists and 
only one of the two additional conditions can be obtained in con- 
junction with equation (30). (See NUMERICAL EXAMPLE AND DISCUSSION.) 

(4) Symmetrical velocity diagram. If V z = t z 2 and 
r 1 = r 2 = r, the symmetrical velocity diagram gives 

V 0 ^i + Vg f 2 = Wr 
or 

l 1 + | 2 -©r 2 (31a) 

Differentiating with respect to r yields: 



it'. 


+ ^ = 20 * 


(31b) 


If r p 4 r 2 1 4 V a , 2 , the symmetrical velocity diagrams 

may be defined by 

t ± (r x ) + i 2 ( r 2) = wr i 2 


(32) 
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Then 

d ^l ^2 

clr-^ + dr-^ “ ^ r i ( 32a) 


3* v e,3 = Y e,l’ ^ v z,3 = y z,2 = Y z,l> the symmetrical velocity dia- 
gram gives very similar flows through the rotor and the stator, 
experiencing the same turning. The change in static pressure or 
enthalpy is the same in passing through the rotor or the stator 
and the stage is therefore often referred to as the "50-percent 
reaction stage." Reference 6 shows that the "blade -profile loss 
is a minimum with the symmetrical velocity diagram if the drag- 
lift ratio is constant. 

(5) Constant axial velocity over the "blade height: 


^z 

3F" 


= o 


(33) 


At very low speed of gas flow with no change in density, the specific 
mass flow is also constant over the "blade height; therefore, there 
is no radial flow across the blade row and equation (26) reduces to 


3r “ r 2 



(34) 


And in case 


<SH 

35 = 


o 


then 



The equivalence of equations (33) and (34) "breaks down, however, 
for current aircraft applications, in which cases the speed of gas 
flow is high. 

If equation (34) is substituted into the radial- equilibrium 
equation (26), there is obtained 
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St z Sr r 

sr - sr - 0 (34a > 

The left-hand side of equation (34a) is the tangential component of 
fluid rotation V X Y; thus equation (34) is a condition for 
potential flow in the free space "between blade rows. 

If it is desired to correct for the effect of the "boundary 
layers at the inner and outer walls of the gas passage, instead of 
equation (33) an appropriate axial- velocity variation close to 
the actual one may "be prescribed in design: 

dV, 

=cp(r) (33a) 

(6) Constant specific mass flow over the blade height. In 
order to avoid radial movement across the blade row in compressible 
flow, it has been suggested (for example, in reference 7) that 
constant axial velocity be replaced by constant specific mass 
flow 


d^ d^. 

dr " dr ~ dr = 


(35) 


Eadial displacement can also be prevented by the use of two condi- 
tions instead of three 


(JG-l dGg dGg 

dr dr “ dr 


(36) 


For designs using either of these two conditions, the simplified- 
radial-equilibrium calculation is more correct. Designs employing 
no radial flow have the advantage that the calculation does not 
involve any radial displacement across the blade row and that the 
two-dimensional-cascade data can be direotly applied. The equa- 
tions for tangential- and axial -velocity distributions derived 
from these equations (equations (35) and (36)), however, are 
difficult to solve and the conditions are incompatible with tapered 
passage in a multistage turbcmachine . 

(7) Constant Mach number relative to blade. For a fixed ratio 
of exit to entrance velocity relative to the blade, the temperature 
or pressure ohange of gas across the blade row at any radius is 
proportional to the square of the entrance Mach number of gas 
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relative to the “blade at that radius. It might therefore be 
desirable to reach the Mach-numb er limitation at all radii. That 
is, for the rotor. 



T Z; i 2 + (urj. - V eA f 


(7-1) 




(37) 


where is the limiting inlet Mach number for the blade. Dif- 
ferentiating with respect to r-^ gives 


(7-1) Mj^ 2 






cor n - 


d ^i 

dr-, 


(37a) 

For stator 





(7-1) 



Differentiating yields 


fr- 1 ’ V S - i 1 + ¥ <k (T */ + Y e/) 


(38) 


(38a) 


Equations (37a) and (38a) may be combined with equation (26) to 
eliminate H. 

(8) Constant turning. If the maximum work for a given size is 
desirable, the limiting turning value may be reached at all radii, 
which will give, of course, a radial gradient in total enthalpy 
after the first blade row. This gradient may not be serious in a 
single- or two-stage unit, but it may not be desirable to use it 
for all stages of an eight- or nine- stage unit. For compressors, 
the limiting turning may be expressed by (reference 2) 
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frjtSJl-Zgii , ( p( g) 
T z,l 

(39) 

or 


Y e ,2 " v e,i +cp(a) v z,i 

(39a) 

where a is the solidity of the blade element at radius r. 
ferentiating yields 

Dif- 

aY e ,2 a*? aT e,i , y ace(2l , . 4T ^1 

dr 2 drp dr^ V z>l drj_ * dr-|_ 

(39b) 

(9) In multistage machines, similar variation in either tan- 
gential velocity, axial velocity, or specific mass flow may be 
specified at the similar stations of each stage 

d£q ^3 3*3 
dr x “ dr 3 dr x 

(40) 

dV z,l dT z,3 ^3 

dr^ dr 3 dr-^ 

(41) 

or 

dG^ dGg di^ 

(42) 


Stages of multistage machines designed for similar variations of 
gas properties from stage to stage are teimed "typical stages." 

Types of design. - A large number of different types of design 
may he obta ine d by different combinations of those conditions 
specified in equations (26) to (42). These designs may be conven- 
iently divided into two groups. In the first group, the condition 
of constant work at all radii is specified In the design. That is, 
equation (28) is specified, which gives: 
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and 


^1 _ ^2 

<^1 ~ ^1 

> 

cLH-j_ 0LH2 dH^ 
dr^ ~ dr £ ~ dr-^ 

J 


(43) 


In cases where the symmetrical velocity diagram is also specified, 
equations (32a) and (43) give 


and 


4 £l ^2 

dTp ~ top “ Wr l 

> 

dH-^ dHg dH^ 

^1 = ^ = ^ ^ 


(44) 


In the second group, the condition of constant work is not specified. 


The following tables present a few types of design in each of 
the two groups. The way in which the degrees of freedom are taken 
up at each station and the known characteristics of each type are 
given. In the tables, subscript i refers to any station in the 
machine; k refers to any station in front of rotor, that is, 
stations 1, 3, 5, „ .; and l refers to any station in front of 

stator, that is, stations 2, 4, 6, . . .. A typical stage may be 
considered as composed of either stations 1, 2, and 3 or stations 2, 
3, and 4. 



aaom x 


*ype 

Sta- 

tion 

Oondlticna 
epeolfied at 
throe ita- 

tlrtna rtf tiny 

stage 

Ono more con- 
dition speo- 
ifiod at oar- 

-Kn-t n jrtotioH 

In anohlnn 

Cfcanaoteristioe of any stage 

1. Jree vortex 
nnl con- 
stant total 
enthalpy 

1 

2 
5 

4tx 

Ccastont work 
dr^ 4r]_ 

dH, 
— 1 

4 *i d? 2 dt» /3v x ar r \ 

5J“EJ’Ei" 0 V " V r i - 0 

Small radial gradient In axial Telooity for 
dH^ flTTo 'IEt rj^rrr* w.nH >i 1 a flow 

35T " “ 3rT “ 0 Constant axial ralooity orar blade height for 

inoooproasibls flow 

2, SjjiarbrloaLl 
▼•looJLty 
diagras and 
OQartOJlt 

total 

enthalpy 

1 

2 

5 

8j™»trloal 

Telooity 

diagram; 
oooartaxrt to* 

ill 4lj 

fTff. 

7^ ‘ 0 

4 ^1 4 ^2 4 ^ ... „ wr l ®t ^l.t ^l.t . 

or x orj ar x -* e,x B * *1 1 v^ei^Tai 

> Tortox-type 

ffi i “a ^ . - /“* 1 e t °i.t r i.t\ r i 1 

a 5^ " JEJ “ 0 v fl,8 - \T + “ fi-*“ ~^) rj, J ^elooltles 

Large nsgatlTe radial gradient <aT axial veloo- 

ity sit s3Jl mTTi vluS3 

5. Symaotrloal 
Telooity 
diagram and 
ccnetant 
axial rol- 
oolty ija 
front of 
rotor 

1 

2 

3 

( Syntotrioal 
Telooity 
i diagram; 

1 ocmatant work 

dT- , dT_ . 


^l d ^2 ... „ “*1 8 t °l,t *l.t 

dr! ’ Txl " “*1 7 e,l ’ T T 2 - 

dEi d^ 4% /tar! 8* rj. r x 

7^ " 7^ * Ef y e,B ‘ ^T* + — 5 — ri ; r 2 

dY £ ! dTj 3 Total enthalpy increeaes from hub to tip 

• ■ — r-- 1 • « 0 Behind rotor, axial Telocity decxoaaaB frc» huh 

to l to tip 

4. DJOB9 tuflijU 1- 
Telocity 
dingroa ht*i 
ooncrtant 
axial Tel- 
ocity la 
front of 
stator 

2 

3 

4 

^z.Z ^1,4 

^ C3 ■ , f i ^ , 

<^2 4r z 

Syi— trloel 
Twlooity 
diagram; 
ocmjb uant vote 

j 

1 


^3 4 ^4 ^ T ttr 3 C t qg.t %t 

^ n ^ n «r 5 7 9,S - ~r - 2^ 

dHg dHj aE 4 _ /tors 8 t l% it r 3<t \ rj 

T e,< ■ \T" + — r 4 

«7, g dT t ^ Total enthalpy inoroaaae from hnh to tip 

■ ij ■ i* f « 0 Behind stator, axial Telooity lnoreaoee froa hnh 

to 2 “4 to tip 
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GROUP I COMCIxmSD 


• i 
i 



I 

I 



Sta- Conditions 
tion specified at 
three sta- 
tiona of any 
stage 



5. Constant 
axial 
Yelooity 


6. Constant 
speolfio 
mss flow 


Senna raria- 
tlcn In 
epedfio 
Bass flow 
and oanstaat 
total 
enthalpy 


Characteristics of any stage 



4tl ^2 

dr^ " dry 

^i-0 “l ®Z 


a7 »,l W t,2 
to l ' ^2 


«V 


dpT -i (lHo dHq> 

3^ -3^° a?£ 

aGy dG2 dflj 

«*1 ^2 “5 


diy dt2 

dH^ dHg dHaj 
dOy 402 dG3 

S^l a lEl a -SFl 


In case of inccttpreaelble flow, this type 
requires no radial flow aoross blade rows 
Tlider sinplified-radial-equillbrltsH oon- 
aldaratlon this type Is equiralent to 
first type In group' 


So radial flow acroee rotor and atator 
blades In case of nontapered passage 
Initial stage is necessary If entering 
gas has different radial gradient in 
total enthalpy than- is required 


Bo radial flow across rotor and stator 
blades in oase of nontapered passage 
Initial guide mne is necessary If 
entering gas has a different radial grad- 
ient cf epeolflo sbbs flow than Is required 
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OJ 

C\5 


®0UP n 


Type 

Sta- 

tion 

Conditions 
specified at 
three sta- 
tions of any 
stage 

One more con- 
dition spec- 
ified at cer- 
tain station 
in maohlne 

Characteristics of any stage 

1 . So radial 
displace- 
ment across 
blade rows 
throughout 
maohlne 

1 

2 

3 

“1 0 
TT— a 0 

dr L 

dOp 

■£r - o 

SB 

35 equal to 

a prescribed 
value at only 
one station 
in the 
maohlne 

<3XJl dQg dQg 

“ 3»b " ^3 “ ° 

No radial f lov aoross all blades in case 
plus 52 equal to a °f naatapered passage 

dr 

prescribed value at 
only one station in 
the machine 

2. Constant 
axial Tel- 
ocity and 
synBetrioal 
yeloolty 
diagram 

i 

1 

1 

2 

3 

fajmsetrical 
Telooity 
J diagram 

Rl «,, 2 

av z i 

B 0 

dT E,l dT z,2 dY z,3 

dr^ “* dr 2 “ dTj " u Coder almpllfied-radial-equilibrium con- 
sideration for.nontapered passage, tan- 
* * gentlal Tolooities and vork done are pro- 

d ^l . “£2 portional to radius (vhpel-type rotation) 

dr^ dr-|_ “ 2ccr]_ and square of radius, respectively, (See 

last part of appendix C.) 

l^x m T*r 

dT ..l w .,s 

dr^ “ dr^ 
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GENERAL METHOD OF CALCULATION 
FOR GIVEN TIRE OF DESIGN 

In. the previous section, examples are given to show how the 
different types of design are obtained by specifying different con- 
ditions compatible with the degrees of freedom available. In the 
calculation of the variation of the gas state through a turbomachine 
of any given type of design, the effect of radial motion has to be 
included, in general, to get a more accurate value. The following 
procedure of calculation is suggested: 

The first step of the calculation is to obtain an approximate 
solution by assuming that the gas flows on cylindrical and conical 
surfaces for straight and tapered passages, respectively, and by 
neglecting the effect of radial motion on the radial- equilibrium 
equation. This method is the usual simplified-radial- equilibrium 
calculation and gives an approximate solution for the case of very 
small blade-row aspect ratio. 

The second step is to obtain the solution for the two extreme 
limi ting cases of zero and infinite blade-row aspect ratio. If the 
difference between the two cases is not large, the solution for the 
case of a certain finite blade-row aspect ratio may be estimated 
from the two extremes. If the difference is found to be large, 
calculation has to be made for the finite aspect ratio. 

In the calculation for the case of finite aspect ratio, the 
method developed using a sinusoidal curve gives exact results if 
such curve is prescribed in the design and 1 b believed to give good 
■approximate results in cases where it is not prescribed in the design. 
Even with this simplification, a series of successive approximations 
is required. 


Simplified-Radial-Equilibrium Calculation 
In this approximation, equation (26) becomes 

_ il dSl d7 z,i 

dr^ “ r 2 cEv[ + "^z,! dr-j_ 

Y r — 0 , 


(26a) 

and 


In the case of nontapered passage, 
equation (26a) is equivalent to 


ri = r 2 = r 3 = r, 
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da ± 

ar7 


\ 2 / 


' 6,1 


dr 


(26b) 


When, equation (26a) is combined with equation (27) 


^1 d ^l 

ri 2 to i 


dY, 


+ Y. 


z,l 


Z A dr-, 


+ 03 


dr. 


2 ^2 d ^2 

I " ^ 1 / " * 2 2 **1 


+ Y, 


dY, 


z,2 dr- 


ZjZ 


i2.fk + T .Esa 

^2 dr^ + v z,3 d^ 


(27b) 


For the case of constant work and nontap ered passage, equation (27b) 
reduces to 

V d7z » 1 V dV Z»2 _1 _ y. d ^l 

T z,l dr = T 2,2 dr + p2 *1^ 8 t ^t dr 


v f^S 1 4 d ^ 

v z,3 dr + ^2 U 


5 _ . d ^r 

3 dr ” dr . 


(27c) 


The radial variation of density and pressure is obtained from 
equations (26a), (A5), and (AS) (appendix A) for the case ^ = 0: 


1 

7 


AA. 

P t,i 


’7 d Pi 

dr^ 


1A. 


dY. 


- Y. 


*2* 


r,i drjL 


(45) 


Pt,i 

P t,l 


n 7-2 _ 

7 P± dri 


'izL 


- Y 


dY. 


r,i dr. 


r A 


(46) 


In these two equations, the last term is very small compared with 
the next to last term and may be neglected. It becomes zero for 
the case of nontapered passage and the resulting equations can be 
more directly derived, as is usually done, by taking the approxi- 
mation involved as the use of 
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_i 7 e,i 2 

Pi c r ± 


for the radial component of equation of motion in place of 


Pl dTi 



- T, 


z,i 



'r,i 


dr 


For a given type of design and given design values of limiting 
Mach number and limiting turning, the variations of gas velocity at 
any inlet station of the stage are known. By inserting the radial 
variation of tangential velocity into equations (45) and (46), the 
radial variations of pressure P 1 (^ 1 ) and. density P 1 (r 1 ) are 

determined. By combining the variation of density with the varia- 
tion of axi al velocity, the variation of specific mass flow G-^r^) 
is obtained. The total mass flow across station 1 is then given 
by 


nr. 


l,t 
*l,h 


2 *r X G i ^l 


The variation of gas properties at the next station (for example, 
station 2) must be such that the continuity equation is satisfied 


p r 2,t r? l,t 

J 270*2 Gg dr 2 =J 2JO--L d^ 

r 2,h r l,h 

By assuming the value of a gas velocity, for example, V z 2 at the 
hub, the velocity and density variations at station 2 can be deter- 
mined in a similar manner, and total mass flow can be obtained. 

The correct distribution of V z 2 > and so forth, which gives the 
correct value of total mass flow, can be obtained in two or three 
trials. 

When the results obtained are substituted into the continuity 
equation (B5) (appendix B), over a portion of the annulus a. certain 
amount of radial displacement across the blade row quite different 
from that assumed in the calculation is obtained. The result 
determined by this simplified-radial- equilibrium consideration is 
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only close to the case of zero aspect ratio of the blade row, which 
will he considered in the next section. This calculation can there- 
fore he used only as a first approximation to the zero -aspect -ratio 
case and also to give a starting value for the calculation of a 
finite aspect ratio. In appendix C, formulas are given in dimen- 
sionless forms for this calculation for two common types of design. 


Zero-Aspect -Eatio Calculations 

The continuity equation in its integral fore (equation (B5)) 
detereines only the function r 2 (rq) and not values of c)V r /dz . 

Two limi t ing cases will now he discussed for which the evaluation 
of the tere c5V r /5z is unnecessary. If the blade row has an 
axial length sufficiently great relative to the blade height (that 
is, if the blade-row aspect ratio is sufficiently small), the tere 
(BT r /5z) i will be negligible in spite of any radial displacement 
required across the blade row. This extreme situation is designated 
the zero-aspect-ratio case and differs from the simplified- radial- 
equilibrium approximation in that the radial displacement across 
the blade row is properly determined and its effect on the state of 
gas is included in the calculation. For a constant-area gas passage 
without taper or with slight taper, this effect is small; conse- 
quently, a successive approximation procedure starting with the 
result of the simplified-radial-equilibrium calculation can be used. 
This procedure may be outlined as follows: 

1. Using the simplified-radial-equilibrium approximation, find 
Y 1 , V Zj2 , t lf 1 2' H l* s 2' Gi, and Gg as functions of rq. 

2. Find r 2 (r-^ by equation (B5), using the value of G 2 (rq) 
obtained from step 1. 

3. Substitute these values of r 2 (rq) into equations (26) 

pmd (27) with the (SV r /dz) tere equal to zero, to get a second 
solution for V z q, V z 2 , £q, | 2 , Hq, and H 2 as functions of 

r r 

4. Bepeat steps 2 and 3 if necessary, using the value of G 2 (r-j^ 
from step 3. 

In step 1 of each process, the same procedure as in the previous 
calc ula tion has to be followed. That is, the variation of gas 
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conditions at station 2 is such that it gives the same mass flow as 
that in station 1. In the case where there is considerable taper 
at the passage walls, it is better to assume rg (r-^) according to 

the taper to start the calculation rather than to use steps 1 and 2. 


Infinite-Aspect-Ratio Calculation 

The other limiting case corresponds to a blade row with axial 
length negligible as compared with the blade height; this case is 
designated the infinite-aspect-ratio case. The negligible axial 
length does not provide space for any appreciable radial displace- 
ment, hence for a nontapered passage rg may be taken as equal to 

v l> or 


Gl = Gg (47) 

and for a tapered passage, a relation similar to the one that fol- 
lows may be used: 


°2 


= Gl 


r l.t 2 ~ r l,h 2 
r 2,t 2 - *l,h 


(47a) 


Because the change in axial length for a very small change in V r 
is also very small, (&T r /Sz) i does not vanish. Although its 
absolute value does not affect the radial motion because of the 
negligible blade-row axial length, the relative value of bV T /dz 
in front of and behind a blade row is needed to determine completely 
the distribution of gas properties at these stations. If the load- 
ing of the blade is relatively symmetrical or is designed to give a 
sinusoidal radial-flcw path, the curvatures of the radial-flow 
path at the two stations are equal in magnitude and opposite in 
sense. Then 



(48) 


In order to combine this relation with equation (26) in a simple 
manner, it may be assumed that 


Y (*r) . y (**) 

\,l\'SrJ 1 - V z,2 X5T/ 


2 


( 49 ) 
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Then combining equation (49) with equation (26) at stations 1 and 2 


av 




+ Y, 


av_ 


dr 


z,2 


dr 


dE 1 

dr 




dr 


^2 dr -jZ dr 


For a typical stage of a given design, equation (47) or (47a) 
together with either equation (48) or equation (49) will completely 
determine the variation of gas properties at the two stations. If 
the variation of gas properties is completely given at one station, 
equation (47) or (47a) alone completely determines the gas state at 
the other station. 


In appendix D, formulas are given in dimensionless forms for 
two common types of design to calculate the variations of gas prop- 
erties for the two preceding limiting cases. The results so 
obtained will give the limits of the variation of the gas properties 
along the blade height. If the difference is large, it is worth- 
while to make the following calculation for a finite blade-row 
aspect ratio. 


a 


Finite-Aspect-Eatio Calculation 

The results previously obtained for a sinusoidal radial-flow 
path is used. Inasmuch as there is no general analytical solution 
of equations (13a), (14a), (24), and (27) that is possible even for 
simple types of design, the method of successive approximation is 
used. 


The procedure of this calculation is somewhat similar to that 
used in the case of zero aspect ratio. In the case where a tapered 
passage or difference in design calls for an r 3 different from 
rp, first estimate the radial positions of gas particles at sta- 
tions 2 and 3, r 2 and rj, respectively, as f met Ions of rp. 
Then calculate the variation of gas properties at stations 1, 2, 
and 3 by equations (13a), (24), and (27), such that the total mass 
flow is the same at the three stations. Using this result, find 
r 2 (rp) and i^(rp) by the continuity equation (B5) and see if 
they check the assumed values. By interpolating rg (rp) and 
r^(i’p) obtained after a few trials at different rp, the value 
obtained by the next try will be close to the correct value. 

For the general case where Hp, Eg, Ip, and l 2 are deter- 
mined by design as f met ions of rp, a rough approximate solution 
gives (see appendix E) 
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r 2 


" r l = 


( r 2 - r l)< 

1 + A 2 


(50) 


This value can be used as a starting value for more exact calcula- 
tion or may be used as the final value for approximate calculation. 
This approximate value of radial displacement across the rotor is 
obtained in the following manner: 

First take as two separate functions 

A e ( r p) the function (r 2 - r^) of r^ satisfying radial- 

equilibrium and total-enthalpy-charge equations (26) 
and (27), respectively, for a given distribution of the 
other variables 

Ac ( r i) the function (r 2 - r^) of r^ satisfying continuity 

equation (14) for a given distribution of other variables 

It is assumed in this method that the radial gradients in Y z and 
f depend primarily on the magnitude of the radial displacement 
(r 2 - r^) and not on its exact distribution. Accordingly A e (r^) 
is set equal to 


A e (r x ) = y e g (r x ) (51) 

where y is the maximum value of A e and g (r-. ) is a plausible 
© 

form for the distribution of A e satisfying the boundary conditions: 

8 < r l,h> = 8 (*I,t) = 0 ^ 

for > (52) 

8’ ( r l,J =0 > S (% n ) = 1 J 

If A c is calculated for a couple of values of A e , it is possible 
to plot y c , the maximum value of A c , against y Q . A fairly good 
approximate solution might be expected to correspond to the point 
yc = ye* This process can be further refined by varying g(r^) from 
the function originally assumed in the direction of the calculated 
function A c /y c . 



40 


HACA TN Ho. 1795 


Distributions of gas properties calculated from rg (rq) given 
by the approximate relation (equation (50)) has been found to agree 
much better with the value obtained by the previously described 
method of successive approximation than that calculated under 
simplified- equilibrium approximations . 


NUMERICAL EXAMPLE AND DISCUSSION 

The method of calculation outlined in the previous sections is 
applied to the typical stages of a compressor employing symmetrical 
velocity diagram and constant total enthalpy, and a compressor and 
a turbine employing free -vortex and constant-total-enthalpy design. 

The calculation is rendered dimensionless by expressing all veloc- 
ities in terms of U.J., total enthalpy in terms of U^ (l) 2 , and r in 
terms of r^. Because the main purpose of calculation is to deter- 
mine th8 magnitude of the oscillatory radial motion and itB effect on 
distribution of gas properties, a nontap ere d constant -area passage is 
used. Heat transfer is assumed to be zero in the calculation and the 
entropy is assumed to be constant at each station. The change of 
entropy across the blades at all radii is taken equal to that obtained 
from the polytropic efficiency assumed at the mean radius. It then 
gives a radial variation of polytropic efficiency decreasing from 
tip to hub. This variation seems to be in the same direction and of 
comparable magnitude with those obtained from experimental data. 

This calculation does not take into account the boundary layers at 
the rotor drum and the outer casing, and is consequently good only 
for the ma i n portion of gas flowing between them. 

In the comparison of different blade-row aspect ratios in 
each design, in addition to the same aerodynamic limitations, the 
same axial velocity at the mean radius is used. The comparison 
between different cases will be slightly different if another basis 
of comparison is used. 

(l) Symmetrical velocity diagram and constant total enthalpy. - 
Because the difference between zero- and infinite-aspect-ratio cases 

is found to be large in this design, calculation Is made for a 
blade-row aspect ratio of 2. The design constants used for all 
cases are: 



HACA TN Ho. 1795 


41 


Hub-tip ratio 0.6 

Limiting Mach number relative to rotor blade 0.8 

? “ Vq 1 

Limiting value of — * — =r— ■ - ■ 3 0.7 

z,l 

Polytropic efficiency at mean radius 0.9 

0.772 


(The last value results from the use of h/^ *= 0.8 in the 

simplif ied-:ra dial-equilibrium calculation, and is used for all 
cases.) The results of the calculation are shown in figure 4. 

Figure 4(a) shows the distribution of specific mass flow in 
front of and behind the rotor for the different cases considered. 

It may be seen that in all cases except the infinite-aspect-ratio 
case, the specific mass flow G-/G^ increases toward the hub faster 
behind the rotor than in front of the rotor; in other words, 
passing through the rotor the gas moves toward the axis of the 
machine. The magnitude of this displacement is obtained from the 
continuity relation and is shown in figure 4(b). In the simplified- 
radial- equilibrium calculation, it is assumed that there is no 
radial motion, but when the distribution of specific mass flow is 
substituted in the continuity equation (B5), quite large radial 
displacement across the blade is obtained. This kind of calcula- 
tion is therefore not a good one. In other calculations, the dis- 
tributions of gas properties are calculated from assumed radial 
displacements that are to be checked with the displacements required 
from the continuity relation with these distributions, and are 
therefore consistent in themselves. The radial displacement to be 
used in the approximate calculation for A = 2 is obtained by the 
approximate formula (equation (50)) and is about 25 percent lower 
than the correct value. 

The variation of axial velocities is presented in figure 4(c), 
which shows that the axial velocities increase toward the hub in 
all cases but at different rates. The high value of axial velocity 
at the hub before the root allows the use of higher turnings 
throughout the blade height without exceeding the limiting value 

(V 2 - WA ! or aC^ at the hub. It also helps to 
give a more unifora 6ach number relative to the rotor blade over 
the blade height. As a result, this type of design gives a higher 
pressure rise and a higher specific mass flow than a free-vortex 
type of design using the same design limitations. In order to 
utilize this advantage fully, the variation of axial velocity 
should be correctly detemined. 
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The calculation of axial velocity based on simplified radial 
equilibrium gives a result close to the zero-aspect-ratio case, 

■which is also true in the distribution of other properties in this 
calculation, because in the case cf zero aspect ratio, the curva- 
ture caused by radial motion is negligible and the difference in 
gas properties caused by the radial displacement is very small in 
the case of a nontapered passage. 

The variation cf tangential velocities is shown in figure 4(d). 
These velocities in different cases vary in a similar manner and 
the difference of magnitude between them is mainly due to the dif- 
ferent value of 6-fc deter min ed by the different values cf 
V z ]_ h /tr t in the various cases. 

Figure 4(e) shows the variation of air angles entering the 
rotor and stator blades. The difference between the simplified- 
radial- equilibrium calculation and the case of aspect ratio of 2 
is significant throughout the whole blade height. The simplified- 
radial- equilibrium calculation gives an error of about -3° at the 
tip of the rotor blade and at the hub cf the stator blades. This 
difference results in an error in the angle of attach; at the design 
point by that amount and the range of operation is also reduced. 

The variation of ftkch number relative to the rotor blades is 
shown in figure 4(f). The simplified- radial-equilibrium calcula- 
tion gives a nearly constant value whereas the more correct calcu- 
lation shows that Mach number actually decreases more than 10 per- 
cent toward the tip for the case of blade-row aspect ratio equal 
to 2. (This variation is only about half of that of a similar 
free-vortex compressor.) 

Figure 4(g) shews the pressure distributions in front of and 
behind the rotor and the pressure rise across the rotor at dif- 
ferent radii. The difference in pressure distributions may explain 
to a certain extent the difference found between measurement and 
simplified- radial-equilibrium calculation. The pressure rise across 
the rotor is fairly uniform in the case of an aspect ratio of 2 
and is a desirable feature. 

The velocity diagrams at three radii for aspect ratios of 0, 

2, and oo are shown in figure 4(h). If this stage is used as the 
first stage of a compressor, the permissible tip rotor speed of 
this design at st anda rd sea- level conditions is equal to 868 and 
826 feet per second for A = 0 and A = 2, respectively. The 
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specific mass flow per unit annulus area corrected to standard sea- 
level conditions is equal to 41.5 and 40.0 pounds per square foot 
per second for A = 0 and A = 2, respectively. 

(2) Free-vortex and constant- total- enthalpy compressor. - The 
design constants used are the same as in the previous calculation. 

In addition, ^0,1 aru ^ ^6,2 are considered to he equal to ~Wq } 2 

and Wq respectively, at the mean radius. In this type of 
design, the simplif ied- radial- equilibrium approximation is equivalent 
to the zero-aspect-ratio case, because axial velocity is also con- 
stant over the blade height due to the constant H and constant i 
with respect to the radius and the negligible curvature effect 
caused by radial motion. That is, the same values of H^_, Hg, 

£l, V z i, and V z £ occur in both cases and the entire cal- 

culation is the same. (See also equation (E6).) 

Because the radial motion involved in this type of design is 
due only to the compressibility of gas, the difference between the 
zero- and infinite-aspect-ratio cases is not large; hence the cal- 
culation for a finite-aspeot-ratio case is not made. 

The distribution of speoific mass flow in front of and behind 
the rotor is presented in figure 5(a). Even in the zero-aspect- 
ratio or simplified-radial-equilibrium case with a constant axial- 
velocity distribution, there is considerable change in density, 
which requires an appreciable amount of outward radial motion to 
obtain the given design conditions behind the rotor. Although the 
amount of this radial motion is small (fig. 5(b)), its effect on 
the variation of gas properties is not entirely negligible. Its 
effect can be seen in the curves of figures 5(c) to 5(g), which 
are somewhat similar to the symmetrical- velocity-diagram and 
constant-total- enthalpy design in nature but of smaller magnitudes. 

If this stage is used as the first stage of a compressor, the 
permissible tip rotor speed at standard sea-level conditions Is 
equal to 758 feet per second for A = 0. This tip speed is about 
13 percent lower than that of the corresponding case of the pre- 
vious design. The specific mass flow corrected to stan da rd sea- 
level conditions is equal to 38.6 pounds per square foot of 
annular area per second for A = 0, which is 7 percent lower than 
that of the corresponding case of the previous design. 
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(3) Free- vortex and constant- total- enthalpy turbine. - The design 
constants used in the calculation are: ^<1, ^t/ a l t = 

T e,i,h/ a i,t = °* 8 > v z,i,m/ a i,t = °* 4 > v e,2 = °> polytypic 

efficiency at mean radius equal to 0.87. For the simplified- radial- 
equilibrium approximation or zero aspect ratio, Y z ^/Ut is constant 

enrl so is Y z "which is found by the continuity relation to be 

equal to 0.87^. The same velocity at station 2 is used for the case 
of infinite aspect ratio, thus making the only difference at sta- 
tion 1. The results of the calculation are shown in figure 6. 

The distribution of specific mass flow in front of and behind 
the rotor is shown in figure 6(a). Because of the constant axial 
exit velocity, the specific mass flow is constant behind the rotor. 

Except for the case of infinite aspect ratio, there is an inward 
radial motion of gas in passing through the rotor (fig. 6(b)), the 
magnitude of which is about two and one-half times that in the 
previous free- vortex compressor (fig. 5(b)). 

The variation of axial velocity in front of the rotor is shown a 

in figure 6(c). An increasing axial velocity toward the hub of 
about 15 percent would be required for an aspect ratio of 2. 

Figure 6(d) shows the radial variation of gas angles entering 
rotor blades. The difference is only important at the hub. In the 
actual case of an aspect ratio of 2, the simplified calculation 
would give an angle of attack 3° to 4° too high at the hub. 

The absolute and relative Mach numbers of gas in front of the 
rotor are shown in figure 6(e). In the actual case of an aspect 
ratio of 2, the Mach number at the hub is about 3 percent higher 
than the simplified calculation. 

Figure 6(f) shows the pressure distribution in front of the 
rotor. For an aspect ratio of 2, the pressures at the tip and at 
the hub are about 2 percent higher and 3 percent lower than the 
simplified calculation, respectively. 

The velocity diagrams at three radii for the zero and Infinite 
aspect ratios are shown- in figure 6(g). 
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SUMMARY OF ANALYSIS AND CALCULATIONS 

In axia l-flow turbomachines, radial motion of gas occurs 
■because of the tapering of the passage and the variation of gas 
conditions across blade rows specified in the design. The direc- 
tion and the magnitude of this radial flow depend on the type of 
design, the tapering of the passage, the hub- tip ratio, the blade- 
row aspect ratio, and the speed of gas flow. Even in the case of 
free- vortex type of design employing a nontapered passage and 
requiring no change in velocity distributions from stage to stage, 
there is an appreciable amount of oscillatory radial motion within 
the stage. 

This radial motion gives an additional texm to the ordinary 
radial- equilibrium equation. In the free space between blade rows, 
this additional teim is very nearly equal to the product of the 
o square of axial velocity and the curvature caused by the radial flow. 
Depending on whether the curvature is positive or negative, the 
radial-pressure gradient caused by the whirling motion of gas is 
decreased or increased, respectively, by this additional term. 

The determination of this radial-flow path requires a long 
process of step-by-step calculation. It is found, however, that 
a sinusoidal radial-flow path gives an effect on the radial varia- 
tion of gas condition between blade rows as small as possible 
without discontinuity in the curvature of the streamline. Inasmuch 
as it represents the major harmonic of the radial-flow path that 
may exist in any design in which the radial blade force is effec- 
tively small, the calculation based on this simple radial-flow 
path gives good approximate results. It probably underestimates 
the effect because of the neglect of higher harmonics. 

The analysis made of the maximum compatible number of the 
degrees of freedom in specifying the radial variations of gas 
properties in stations between successive blade rows of a turbo- 
machine shows that under the conventional design procedure the 
designer is free to specify two conditions at any one station and 
one condition at each of the remaining stations. The various ways 
to use up these degrees of freedom and the resultant types of 
design obtained are discussed. 

The usual method of calculation, which neglects the radial 
motion, gives results close only to the case in which the axial 
length of the blade row is much larger than its radial length, 
and is not good for the case of a finite blade-row aspect ratio. 
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The difference "between the results obtained "by the usual method and 
the method suggested herein is found to he quite large in a design 
employing constant total enthalpy and symmetrical Telocity diagram. 
Calculation made for this type of compressor using the same limiting 
Mach number, same limiting turning, same axial Telocity at the mean 
radius, and for a blade-row aspect ratio of 2 gives the following 
differences between the usual and the suggested method: 

1. The radial variation of axial Telocity in front of the 
rotor is 13 percent for the usual method and 28 percent for the 
suggested method, and the radial variation of axial velocity behind 
the rotor is 53 percent for the former and 40 percent for the latter 
(all expressed in terms of their values at the mean radius) . 

2. The air angles differ from 1° to 3° at the hub and at the 

tip. 


3. The radial variation of Mach number relative to the rotor 
blade in the former is 9 percent lower than that in the latter. 

4. The radial variation in static-pressure rise across the 
rotor is 13 percent for the former and only 2 percent for the 
latter. 


5. The mass flow in the former is 4 percent higher than that 
in the latter. 

6. The permissible rotor speed in the former is 5 percent 
higher than that in the latter. ' 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 14, 1948. 
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APPENDIX A 

DERIVATION OF EQUATIONS IN ANALYSIS 
From equation (3), the definition of h, and the relations 


and 


there is obtained 


du = c v dT 




From equations (1) and (A4) 

“ = £ d (p) + 4 (ir) 

From equation (2) and the definition of h 


Using equation (A4) 


T ds = dh - ^2 
P 


Tds=-2^ r d(2')-^£ 

7-1 V P / P 


(Al) 

(A2) 

(A3) 


(A4) 


(A5) 


(A6) 


(A7) 


From equations (2) and (3) 
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d /£\ _ _2L dp 

\B/ 7-1 P 7-1 P 


= — ^=- d f^lOgg -2- 

7-1 V P7 


A; d (lose pj- d d°Se p) 


From equations (4) and the following relations 

i-H + < ? - 7 > ? 


and 


- V X (V X T) + |v7 2 = F - ~ VP + £ V 2 ? + ^ V (V -7) 


(A8) 


(A9) 


= T i + ^ + Y *(^ + I 


(A10) 


(V*v) V = -| VV 2 - T X (v X 7) 

there is obtained 

57 
St 

From equations (A5) and. (A7) 

|w2 + i7p = VB-^_7(£) + iv E 

= VH - T Vs 

Combination with equation (A10) yields 

VH = F + TVs + V X (V X 7) - ||.+ £ jj 2 ? + | V(V«V)J . (7) 

From equations (l) and (A6) 

M = T , 1 £E + v.H 

Dt x Dt + PDt + v Dt 
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By combining with equation (4) 


m 

Dt 


= t 55 + P^ +? -{ f + p [’^ + I v(vt)]} (ah) 


When equation (All) is combined with equations (2) and (5) 


m 

ut 


=:Q + p'*^lt + 7 *{ F + p [ v2? + f V(V * Y) } 


From equations (2) and (5) 


T — = Q + ^ 
Dt y P 


or 


M - St . E 

Dt ~ T p 


k 


( 8 ) 


( 9 ) 


From equations (6a), (3), and (A9), the following form of continuity 
relation is obtained: 


V ‘ 7 + 7-1 Dt log e T ~ Dt ( e) = 0 ( 10 ) 

Equation (8a) is obtained by applying motion and energy equa- 
tions to a mass system with a fixed control surface as shown by the 
solid lines in figure 3. Under steady axially symmetrical flow, 
the mass inflow da^ in time dt is equal to the mass outflow dm^ 
in time dt, the state of gas within the control surface is not 
changed, and the state of gas at stations i and j is constant 
with respect to 0. By equation (A10), the sum of the tangential 
blade force and the tangential viscous gas forces exerted by the 
surrounding gas particles on the system is equal to 
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The torque about the z-axis exerted "by these tangential forces on 

D(r Y e ) 

the system is therefore simply — ^ dm^ and the work input to 

the system by these tangential forces in time dt is equal to 


D(rT 0 ) 

Dt 


dm^ codt 



(rT e ) ± ] dmi 


In passing from station i to station i, beside receiving this 
work input, the gas particles are doing work against the axial and 
radial viscous forces exerted by the surrounding gas particles. 
However, this negative work is usually small, and if it is assumed 
that the heat generated from this frictional work is added back to 
the gas stream then the heat addition cancels the negative work, 
and the energy equation for steady flow gives 



(rY e )i]dmi 


or 


del 

Dt 


Q + co 


D(rV e ) 

Dt 


(8a) 


where Q denotes the rate per unit mass at which the gas stream 
sheet is receiving heat from external source through blades or 
other passage walls. 

It may be noted that for steady nonvlscous flow, equation (8a) 
can be obtained by using equation (11).' For such flow, equation (8) 
becomes 


s-«** 


From equation_(ll) it can be 
However, S’ • U is equal to 
, , ’ D(rY 0 ) 

tion (7h), equal to — — . 


seen that F • Y is equal to F • U. 
FgU or F e ru) and F@r is, by equa- 

Hence the preceding equation becomes 
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m 

Dt 


D(rTg) 
Q +a) Dt 


(8a) 


When equation (12) is given, the entropy change can he obtained 
in the following manner. From equation (A9) and (3) 


Dt (r) “ 7-1 Dt ( log e " Dt ^ 1 ° S © 


But hy equations (12) and (3) 


5S lo Se P - ^ lo Se T 


Tvf* lo Se ^ 


(A12) 


Substitution into equation (A12) gives 

£B _ / J. l\ 

Dt “ " l 7-1 n-1 J Dt 

<A13) 

For steady axially symmetrical flow, equation (A13) reduces to 


§ - E ( T r S ^ I + V 2 |r log. t) (13) 

For successive axial stations 
apart, equation (A13) gives 

B J (r 3> - s i (r i) - E 


_ T> a- 7 
(n-1) (7-I) 

Inasmuch as the temperature change between the two successive sta- 
tions is small, the temperature ratio can be considered equal to 
the enthalpy ratio: 


i and J a short distance 


[log e Tj (rj) - log e T ± (r ± ) 


log. 




(A14) 
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_Ll 


% *i 


- Z±! 

% - — 


Substitution of equation (A15) into equation (A14) gives 


(A15) 


Sj (rj) - SjL (r ± ) = B loSe 


rr IsL 

H J ' 2 


% --r 


(13a) 


The density ratio between the two stations is' obtained from 
equations (A9) and (3) 



Combination with equation (A15) yields 



(A16) 


(A17) 


Substitution into equation (14) gives 



( t, 2 VI -fl 

( H J - 2 / e E r 3 ^3 - Y z,l 




1 

i7-l 


e 


fi 

E r i ^i 


(14a) 
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APPENDIX B 


DETERMINATION OF RADIAL DISPLACEMENT 
FROM CONTINUITT EQUATION 

Equation (14) may be written as a linear differential equation 
for r g2 as a function of r^, provided Gg(r^) is known 


d < r 2 ) = ? Gl ^l* 

<3^1 Gg ( r]L > 1 


%t 


r 2 2 = 


°1 (r l } 2 

G^J 2r l *1 + r 2,t 


(Bl)' 


when divided by r g ^ 


r 2 > 

r 2,ty 







°2 (^l) %t r l,t 


If G^ and Gg are not known, and only G i/G^ t and Gg/Gg ^ 
are known, a modification is necessary " * 


“, ty 


1 - 



G! 

^t fri) r l> ' ^1 


^ (X!) 

G 2,t x 


(B3) 

The value of G-^tA^t is found by the condition that total mass 
flow at stations 1 and 2 is the same 
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APPENDIX C 


SIMPLIFIED-RAI1J1AL-BQIJILIBRIUM EQUATIONS 

Equations to calculate distributions of gas properties under 
the simplified-radial-equilibrium approximation for a few types of 
design are given. 


Free Vortex and Constant Total Enthalpy 
For this design 


dIL|_ 

dr7 


0 , 


dr^ 


(Cl) 


from equation (26a) in the section entitled " Simplif ied-Badial- 
Equilibrium Calculation," 


dV 




d r-t 


(C2) 


The variation in tangential velocity is, by equation (Cl), 

r. 


v e,i ~ v e,i,t 


• j,t 


(C3) 


At each station, by using equations (46) and (C3), 

Ft ,.7-2 dp _ (^0,t r t) 

V Sr" r 3 

When the preceding equation is integrated from r to r t , and 
the relation 


a = 



is used, there is obtained 
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This equation holds for all stations, provided the appropriate 
values of (V & ^/a^) are used. It follows from equation (C2) 

that at each station 


-G m £- 

Gt Pt 


(C5) 


The radial position of gas at station 2 or 3 is then obtained 
by n u merica ll y integrating equation (B5) using distributions of 
specific mass flow given by equation (C5). An alternate method is 
to expand the right-hand side of equation (C4) into a binominal 

series . Because — ^ 

0.15, three terms will be sufficient. Let p r represent the 
average density in the annulus between r and r^, then 


fe>) Pf • 


is usually less than 


fr 

Pt 



2rtpr dr 


jt(r t 2 - r 2 ) p t 





as 


r 2,t 


•' 1:'.^ 1 V_ , dr, = 


J r,l v z,l Ui 'l = / 2ltrp r,2 V z,2 ^2 

' <2 


(C6) 
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APPENDIX C 


SIMPLIFIED-BADIAL-EQUILIHKEUM EQUATIONS 

Equations to calculate distributions of gas properties under 
the simplified-radial-equilibrium approximation for a few types of 
design are given. 


Free Vortex and Constant Total Enthalpy 


For this design 


S. ft fii 

dr^ * dr^ 


0 


(Cl) 


from equation (26a) in the section entitled "Simplified-Eadial- 
Equilibrium Calculation,” 


dV 




dr-t 


(C2) 


The variation in tangential velocity is, by equation (Cl), 

V -V 

v e,i " v e,i,t r ± 


(C3) 


At each station, by using equations (46) and (C3), 

„ Pt n 7-2 dp _ ( v 0,t r t) 

^7 3r ' 

When the preceding equation is integrated from r to r t , and 
the relation 


a = 



is used, there is obtained 
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This equation holds for all stations, provided the appropriate 
values of (V 0 ^/a^) are used. It follows from equation (C2) 

that at each station 


1 °. = -£- 

G t Pt 


(C5) 


The radial position of gas at station 2 or 3 is then obtained 
by numerically integratin g equation (B5 ) using distributions of 
specific mass flow given by equation (C5). An alternate method is 
to expand the right-hand side of equation (C4) into a binominal 

series . Because ~- 

0.15, three terms will be sufficient. Let p r represent the 
average density in the annulus between r and r^., then 


(%■) &)'-■ 


is usually less than 



and inasmuch as 
r l,t 


r 2,t 


2ltr ?r,l V z,l ^1 = 


2irr P r ,2 V z,2 ^2 
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or 


p r,l \,1 ^ r l,t 2 " r l 2 ^ = p r ,2 V z,2 ( r 2,t 2 - r ? 2 ) 


’,1 (r l,t 2 ~ r l 2) JzjZ _ p h,l (r l, 

2 _ „ 2\ V_ , ~ 2 


2 

- r 2 ) 
t r l,h 3 


p r,2 ( r 2,t " r 2 ) zA Ph,2 (%t 2 " r 2,h 2 ^ 


hence 


r 2 \ 2 


~2,t/ 


= 1 - 


1 - 

(%hf 

\ r 2 , t / 

y[ 

f *iA* 

\ T l,t) _ 


r l\ 2 




■ -gfa>i 

p r,2 

' Ph,2 


Then 


The change of total enthalpy across the rotor is 

Eg - % = 0)U 2 - | x ) 

h. - 1 „ -h) 


H- 


Hn 


For compressor 


2,h 


*2 , r l,h 7 z,l,h r l,h Ve > 2 > h Ye > 1 > h U 1>t 2 

TT TT rr I r 

*1 


r l,t u l,t 


'z,l,h 


X 


where 


the quantity 



V 0,2,h “ 
V z,l,h 


f 0A,* 


is to he assigned hy the 


designer. Inasmuch as 
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where is the limiting Mach number to be chosen by the designer; 



(C8) 


The pressure distribution at each station is obtained by rais- 
ing its density distribution (equation (C4)) to the power 7. The 
pressure changes between the stations at different radii are 
obtained by combining these pressure distributions with the pres- 
sure change across the rotor at the radius where the value of the 
polytropic exponent is known or assumed. The angle that the gas 
velocity makes with the axis of the machine at any radius is 
obtained from the known tangential and axial velocities. 
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Symmetrical Velocity Diagram and Constant Total Ent halp y 


For the case of nontapered passage, r^ = = r and from 

equations (28) and (32) 


dr ~ dr 


cor 


-\ 


f - ^ r l»t S t u l,t 
5 1 2 2 


, cor 2 r l,t ^t fl.,t 
f = + — ^ 5 


> 


(C9) 


or 


ILi 

1 

r 8 t 

r l,t 

u l,t 

2 

r l,t 2 

r 

v e,2 

1 

r & t 

r l,t 

U x,t 

2 

r i,t * 2 

r 




> 


(CIO) 


By substituting equation (CIO) into equation (46), 


y ft 7-2 S£ _ cofr _ 5 t P l,t 2 1 

P-j.^ Sr ^ ^r ^ 




where the minus sign is used for station 1 and the plus sign for 
station 2. Integration from r to r t yields 


7 

7-1 



.7-li 



r 2 ) 


5 t U l,t 


lo g e 



r 


1 

8 


< r i,t »t u i,t> 2 




or 


C 
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7-1 °t 2 



1 



rfc 


45 t 




(Cll) 


where the plus sign is used for station 1 and- the minus sign for 
station 2 . 


The variation of axial velocity is obtained from equation (26a) : 


'z dr 



i 


*1 

dr 



r l,t 5 t U l,t\ 
2r ) 



(C12) 


Integration from r^ to r 


gives 







(C13) 


where the plus sign in the last term is used for station 1 and the 
minus sign for station 2. Equation (C12) also gives 


dV, 


z,2 ar 


2 il 


- Y, 


<3Y. 


dp 


Z A 


5 t u i,t 2 


(C14) 


In the case of a tapered passage, the gas is assumed to flow 
in conical surfaces, which gives the value of r 2 as a function ' 

of r-j_. Equations (28) and (32) give the distributions of | and 
Y q as shown by equations (D6) and (D7), respectively, given in 
appendix D. The distributions of axial velocity and density at 
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station 2 are the same as those given by equations (D8) and (D5). 
After these distributions are known, the distribution of specific 
mass flow G is known and the radial displacement is found from 
equation (B5) . 


In compressors of this design, the maximum value of 


V e 2 - V 0 i 
r l,h g * 2 9,1 

Y *,l 

set by the designer. Then 


is usually at the hub. Its value there is to be 


and 


8 h = 


t 


2-i_b 


I 


r l,h u l,t 


l^h = r l,h 


2,h _ _ v 

V e,2,h v 0 ,l,h v. 


z,l,h 


z»l.b 

D i,t 


2,h 


5 + = 
t , r. 


l,h g, , l,h z,l,h r l,h 

* x i,t h r l,t U l,t ? z,l,h 


6,2, h ~ v 6,l,h 


(C15) 


In this type of design, the limiting Mach number is usually at 
the hub. Hence the denominator of the last term of equation (C8) 
should be replaced by 


1 

(7-DMl 2 



, U l,t ) 




The rest of the calculation is the same as in the previous design. 


Symmetrical Velocity Diagram and Condition on Axial Velocity 


For constant work over the blade height and nontapered passage, 
equations (C9) to (Cll) still hold. If equals 0, equa- 

tion (27C) tinder Simplified Radial-Equilibrium Calculation gives 


dV, 


%1 




dr 


= ~2 r l,t 
r 


U l,t wr = 


5 t U l,t 
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Then 


^ - 4 ♦ los e -f 

=A,t v t 1)t 2 e v 


(C16) 


The enthalpy variation for this case is obtained by using equa- 
tions (26a), (27a), and (C9) 


dH x _ dHg _ i &t 2 _ i / r 2 r l,t 5 t U 1,1 

dr dr jJ2 ^ dr r 2\ a) 2 + 2 


cor 


U 2 


r \ 2 . 6 t 


■l,t 


*t - H _1 

tj 2 4 

U l,t 


1 - 


l,t 


+ | & t log e T 


On the other hand if 


dr 


(C17) 


ay 




u. 


'z,2 d.r 


1^ 


2 5, 


v z,2 = 
^z,2,t 


2o u 2 \V 2 

= II - -^4- log. £•) 


V 2 
v z,2,t 


3 e r+ 


(CIS) 


Et - E x 

2 4 


U- 


l,t 


1 - 


*)’ 


- I 5 t lo Se ~ 


(C19) 


In order to make the distribution of E independent of 5^., the 
compromise condition 


av z dV z 2 

T *,l -1^ - - 0 
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may be used. When this relation is combined with equation (27c), 
there is obtained 


- (i + iog e ~ 

V z,l,t V ^Z,l,t 2 t 


1/2 


(C20) 


z,2 


. (i . lt_%t 


z,2,t 


'z,2,t' 


3 -Se * 


1/2 


(C21) 


When equation (26a) is applied to stations 1 and 2, and the two 
resulting equations are added 


Then 



(C22) 


If the requirement of constant work over the blade height is 
abandoned, it is possible to require no change in the axial - 
velocity distribution across the blade rows, which is the incom- 
pressible form of the condition of no radial displacement 


dV, 


z ,l 


dV, 


z,2 


dV. 


dr. 


dr c 


z ,3 


dr. 


= 0 


(C23) 


When equation (C23) is combined with equation (27b) for the case of 
nontapered passage 


ii = li ^2 / d| 2 d|j\ 

r 2 dr r 2 dr \ dr dr J 


(C24) 


By substituting the symmetry conditions (equations (31a) and (31b)) 
into equation (C24) 
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if 

r2 51 dr 


= — (cor 2 - 1 2 ) ^CUr 


- col 2 COr - 2 


cL£; 

dr 


dr j 


Simplification yields 


% 

dr 


2 ii 

r 


6l 


r l,t. 


(C25) 


or 


'6,1 _ r 


V 0,l,t r l,t 


From equations (32) and (C25) 

h 


h,t V r l,t 


(25a) 


(C26) 


v e,2 _ r 2 , t ( r 


'e,2,t r i/t V %ty 


(C26a) 


Thus wheel-type rotation exists both in front of and behind the 
rotor. By equation (27), 

ABg ^1 _ / d ^2 _ ^l\ 

dr dr “^ydr dr/ 


2cor 


r l,t 


2 ^2,t " ll,t) 


- 2fl*Jl,t 5 t 



(C27) 


which means that the work done is proportional to the square of the 
radius . 
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APPENDIX D 

METHOD OE CALCULATION FOE ZERO AND INFINITE ASPECT RATIOS 

The method of calculations is given for the two types of design 
used in the numerical examples. 


Free Vortex and Constant Total Enthalpy 

Zero aspect ratio. - In this type of design., the zero-aspect- 
ratio case is the same as that of the simplified-radial-equilibrium 
approximation . 

Infinite aspect ratio. - By equations (47) and (A17) 




- < 7 e,x £ * 7 *,i 2) 

2% - < T e,2 2 + T z/> 


(Dl) 


An additional relation between V z ^ and V z ^ 2 is necessary 

in order to solve the equation. In the section entitled "Infinite- 
Aspect-Ratio Calculation" two equations- are suggested. For this 
design, equation (48) gives 


or 


*Vl dV z,2 

dri dr 2 


= 0 


V z 1 + V z 2 = constant 


(D2) 
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and equation (49)- gives 


dV, 


z,l dr 




dV, 


+ V, 


z,2 dr 


= 0 


or 

+ ^•z.,2 ~ cons b an 'k (D3) 

Also from equations (E4) and (E5), neglecting the square term in A e 


7 zl V Z) 2 ~ content (D4) 

The three preceding equations give practically the same results. 


A convenient procedure of calculation is as follows: 

(1) In order to compare the result with other cases, the same 

value of v z d m U8e(3 -* F rom equation (Dl), V Zf 2,m iB 

determined. 

(2) Insert these values in either equation (D2), (D3), or (D4) 
to ohtain the constant in the equation. 


(3) Assume a number of values of Y z ohtain Y z 2 "by the 
same equation. Then use the following equation, which is obtained 
from equations (Dl) and (C3), to solve for r/r t 




( v e,2,t) 2 






(4) Plot Y Z)1 
values at the values 


and V Z) 2 
of r/r t 


against 

desired. 


*/*t> 


and obtain their 
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When the distribution of axial velocity is known, the density 
variation at any station is obtained by applying equation (A17) at 
the station 



(D5) 


The pressure variation at each station is obtained by raising equa- 
tion (D5) to the 7 power. The pressure changes across the stage 
at different radii and the air angles are obtained in the same 
manner as in the simplified-radial-equilibrium calculation. 


Symmetrical Velocity Diagram Plus Constant Total Enthalpy 

Zero aspect ratio. - With radial displacement not equal to 
zero, the equations for tangential velocities are different from 
the expressions of equation (CIO). From equations (28) and (32) 


d£, d| 2 

-r- : - = ~ - = car-, 

drq dr-L 1 


cor 2 r, A 5, U. 


i = _i- - 
S 1 


l.t t l,t 




® P 1 *l.t 5 t U l.t 
+ — 1 


2 2 


(D6) 


and 


Y e,i _ ii _ r i _ r i,t 5 t 

u i,t r i u i,t 2 r i,t 2r i 




' 0,2 _ 


r i,t 8 t\ r i 


u 


l,t " r 2 u l,t " \ 2 r l,t 2r l / 


(D7) 
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From equation (26) -with the last term negligible for this case, 
inasmuch as r 2 - rq is much less than L, 


d7 z.i = _ -i- t ili 

I 1 drq r± 2 i arq 


For station 1, from equation (D6) 

2 


OT Z .l 1 r 1 r l.t St Mm- tt 2 (l r l S t 

,1-Ef-^TS * * 2ri 


Integration from r^ to rq yields 


’Vi^ /v - 

U 


i,t/ V u i,t 


+8 t lo Se— W ^ 


?h2L 


V r i,t r i,t 


which is the same as equation (C13). For station 2, 


dV. 


z^2 


'z,2 dr. 


^z,2 *^z,2 




r 2 

r 


2r. 


\ (0)2r l 2 + 6 t U l,t 2 ) to i 


2 n j * 




'l,t. 


l,ty 


Integration from rq ^ to rq yields 

r l 


%2j . /Iz ^ hj 


r r l,t 


^h 

r l,t 


fif /IfL + A-W - 1 

v l r i,t r i > lr 




r l,ty 


which differs from equation (C13). 


(D8) 
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The density distribution at station 1 is the same as the 
simplified radial approximation, whereas that at station 2 is 
obtained by using equation (D5) . The solution of this case is a 
process of successive approximations. Yalues of rgCr^) obtained 
in the simplified radial approximation can be used here as the 
starting values. Then the distributions of Y q 2 > V z 2 > 
and Go are 'calculated from the preceding equations, and new values 
of r 2 (r 1 ) are computed from equation (B5). Usually, only two or 
three cycles are necessary to obtain the correct value, as the 
difference between this case and the simplified-radial-equilibrium 
approximation of this type of design is small. 


Infinite aspect ratio. - The first equation for the condition 
Gq = G 2 is the same as equation (Dl). The second equation 
necessary to solve this case is a little more complicated than that 

in the previous type of design because ^ ^ 0. If equation (49a) 
is used 


dY z q dV_ o 

z,l dr z,2 dr 


d| x 



= - ^ (lq + tz) ° “d) 2 r 

Then 

Y 2 + Y„ o 2 = - 0> 2 r 2 + constant 

z,± Z,£ 

or 

Y 2 Y 2 2 

~- L 0 — b - - Z - ? - 2 - - = - + constant (DIO) 

U t 2 U t 2 r t 2 

In order to compare the result of this case with other cases, the 
same value of Y z ^ m may be used. Then from equation (Dl), 

Y z 2 m 'l® f oun( l> and the constant ' in equation (DIO) is evaluated 
by using this set of Y z ^ m and Y z ^2,m' A f®** values of Y z ^q 
are assumed at any other given radius, with corresponding values 
of Y z 2 obtained from equation (DIO) . The correct values of 
Y z q and Y Zj 2 that will satisfy equation (Dl) are obtained by 
interpolation . 
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After the distribution of axial Telocity is known, the density 
distributions are obtained from equations (Cll) and (D5), and the 
pressure distributions, total enthalpy change, and air angles are 
obtained in the same manner as before. 
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APPENDIX E 

APPROXIMATE VALUE OF RADIAL DISPLACEMENT ACROSS BLADE ROW 


FOR GENERAL CASE IN WHICH H x , Hg, l 1} AND % 2 
ARE DETERMINED BY DESIST AS FUNCTIONS OF r x 


In the latter part of appendix D, distribution of axial veloc- 
ity is expressed in terms of known H, £, r^, and r 2 (r]_). 

Alternatively, this distribution can be expressed in terms of radial 
displacement and its value determined by the simplified-radial - 
equilibrium calculation, for which A e = 0. For the case of 
nontapered passage, it is seen from equation (21) 


dV, 


z,l dr 




dV, 


= V, 


z,l,s dr 


z ,l,s 


(f) ! 


Z,1 


(El) 


and 


■ ? -V „ dYz ;. 2 .u.?._^ (if f v 

z > 2 dr-j^ z,2,s 2 \JL/ \ ■ z t 2 J dr-j^ 



By substituting (r-^ + A e ) for r 2 , expanding 

binomial series, and neglecting terms of greater 
equation (E2) becomes 



order than 


(E2) 


in a 



dV, 


z,2 dr n 


z ,2 


dV. 


= V. 


z,2,s dr-. 


z,2,s 



+ 


. a fe 

^2 d^ 





(E3) 



-0 

CO 


If a 0 (r]_) 1 b known, equations (El) and (B3) may "be solyed as linear first-order differential 
equations in T, ^ and V„ ? 2 , respectively, giving (omitting the subscript 1 on r) 

“i— "i" ' 


■ r i 



and subscript' m may here refer to any radius between hub and tip. 


- I 

- 1 


i 
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For the limiting case of zero aspect ratio, the last term in 
equation (El) approaches zero, so Y z -j_ = Y z i^ 8 , whereas Y z 2 
is obtained by integrating equation (E3) with the third term 
neglected. 



(V. 


z,2,s,m 


- V, 


z,2,m 


•) + 2 




dr 



(E6) 


When equation (E4) is integrated by parts and A e (r-^) is 
replaced by y e g(r) as in equation (51), there is obtained 


▼z,l = 7 Z,1,S + 0 


(l) *1 ( r) 


(y , 2 - V n 2 \ 

^ v Z,l,m v z,l,m,s J 


- e 


(|) ye^fr) 


r r 


V z,l,s dr 0 


-{if *1 M 


dr 


■■m 


in which 




Vj. (r) = 


g(r) dr 


‘•m 
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If it is desired to compare the general case with other cases on the 
basis of the same V z ^ m , then V Z}1}iaL = V Zf l,m,B a By the use of 
the mean value theorem of integral calculus, the preceding equation 
can he written as 


'z,l “ V Z,1,S 


= - e 


(l) y e ^ 


= V, 


z,l,s 


z,l,s 

{if ^ (r) , 

0 - X 


-(f) 2 Je *1 < r > 

5 - 1 


(E7) 


m> 


where T_ „ is a mean value of V- t s between r and r^ 

ZtyJ-jD " j w 

the mean depending on the choice of the function g (r) . If the 
approximation is made in letting V z , 2 ,m = ^’z^2,m,s^ equation (15) 
may be written as 


rr 2 v 2 _ _ 2 

*z,2 “ y z,2,s “ v z,2,s 


-(f) 2 y e <*>2 m 

j ' - 1 


+ y e ^3 ( r ) 


- y/ cp 4 (r) 


(18) 


where 


(r) = 


= «P X (r) + i y ( 


<P 3 (r) = e 


y e 

-{if * (r > 


1 

2 v e 


[1 

r. <r)i 2 - r 

-i 2 

S (r m )J k 

L 

- 


,2 ^\2 
rr 


(f) * < r > 




m 


and 


<P 4 (r) 


-dfcpM p (f) 2 cp (D ^ ]2 

- 6 1 e 4 dr ^2 


dr 


-'r, 


m 
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The change In the distributions of Y Z)1 and Y Z)2 vith the im-rlmum displacement y 0 for a 
given g (r) is now d etermined by differentiating equations (E7) and (E8), assuming that r m , 


V_ i . Anti V ~ nrv? 1nflflnflnflan+ nf it 

7 Z,C t B “ •'e 




2 4 T, (1 - (if ^ to 


(f) ^ 'Pi to 


(E9) 


2 loft, V, 9 - 

“Je ~ 


79 ” 2 (r Y*, 


h, 2 2 W 


V‘ 


(r) + y„ < 


g(r) 

* 


*’ m u i i 


v 


y e *P3(r) cp 4 (r) • y 0 2 dcp 4 (r) 


o + o — — — - 2y, 

V ^ v ^ 2 v Zf2 


o j 

* “jr 


(E10) 


By subtracting equation (E10) from (E9) and neglecting three small terms containing y e /V e 2 


d V z,l 1 /k\ 2 m . J ’z.l.a i 

dy Q l0 Se Y Z}2 ~ 2 (l) ^ (r ^ y' ,fc 0 


(e) y e ^ GO 


'ajl 


. _g»2»a 
V 2 

T(7 O 


•(e) y 8 ^ to 


1 + 7- 


[s(r )] 2 - [s(r m )] 2 

cp-, (r) 


i 

/ 


1 *5 GO , , 

£ 2 (E11 ) 
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The equation of continuity, equation (14), may he written as 


PjVz,! 
p 2 7 z,2 = 


, y c s( r )’ 

~ — 
d y x 

i + 

r 

1 + 7 0 S M 


(E12) 


by replacing r 2 with. r^_ + y c g (r). It is here "being assumed 
that the displacement A c for the continuity equation has the same 
form as, "but may differ in magnitude from, A 0 . If the variation 
of P}/p 2 (r) with y c is neglected, that is, the density distribu- 
tion is assumed to be determined primarily by the tangential veloc- 
ity distribution, differentiation with respect to y c for a 
given g (r) gives 




(E13) 


If the same distribution of V z 1 and Y z 2 satisfies both the 
continuity and equilibrium equations, y c is a function of y e 
determined by the differential equation, which is obtained by 
dividing equation (Ell) by equation (E13) : 



<1 


+ y e 


[s(r)] 2 - [g(r m )] 2 

cp^r) 



(E14) 
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la order to evaluate equation (E14), the form of g(r), which is 
implicit in the equations used, must be found. An order- of -magnitude 
result may "be obtained, however, by equating the right-hand side of 
equation (E14) to a constant -K^ and determining the value of 
the constant from the boundary conditions on g(r) . Because this 
assumption involves setting dy c /dy e equal to a constant, it is 
equivalent to the assumption already stated that for the selected 
A 0 = y 0 g(r), the corresponding A c differs only in amplitude. 

In order to obtain the order-of -magnitude result, the right-hand 
side of equation (E14) is simplified by 


(a) Setting the first two terms in the bracket equal to 2 

(b) Considering the terms involving y c and y Q negligible 

when they are compared with unity 


(c) Ignoring the last term in the bracket because it contains 
V z ,2 2 in ill© denominator. (If equation (E14) is 


written in terms of r/r^ instead of r, the term 
\2 /jtr-fcN 2 

-=— I , which is about 250 for A = 2 


(!)’ 


—J becomes 
and ~ = 0.6. ) 


L J 


As a result of this simplification, equation (E14) becomes 


d ?c = (f) ’Pi < r > ' ^ 

** = S & + £ g(r) = 


Rewriting equation (E15) gives 


(E15) 


b 



When equation (E16) is differentiated with respect to r 
relation (r) = g(r) is used. 


s(n) + ~ s(r) + 


dr 




0 


(E16) 
and the 
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This equation gives g(r) as a Bessel function of the first order 
and argument (jt/lK). The value of (n/LK). and thus of E, is 
determined h y the "boundary conditions = g(r^) =0. In order 

that g(r) have a single maximum, the first eigenvalue of this 
boundary-value problem must be taken. A satisfactory approximation 
to this solution may be obtained without involving Bessel functions by 
replacing g(r)/r in equation (E16) by g(r)/r m , differentiating, 
and solving: 

• 

r 

“ 2r 

g(r) = e m ( % cos Zr + K2 sin 2r ) 


where 


Z = 



1 

(^m ) 2 


The boundary conditions determine, using the first eigenvalue 
for K, 


K 


ji 

*V r h 


2r™ r-r>, 

g(r) = K 3 e m sin * (E17) 

0 r t~ r h 

and therefore, 

(E'f JL = * 2 1 ?t 2 

UJ K 2 (r t -rn) 2 4r m 2 (r t -rh) 2 


(This approximate equality is correct within 1 percent for 
(r^/r^) > 0.5.) Substituting this result in equation (E15) gives 


^e 


- K 2 



(E18) 


In this very rough approximation dy c /dy Q is therefore equal to 
minus the square of the aspect ratio. By integrating equation (E18) 


\ 
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and "by letting y c s equal the value of y c corresponding to 
y e = 0 (simplified-radial-equilibrium approximation): 

y c = y 0 ,a - a2 y e (E19) 

A solution corresponds to y c = y e = y, which when substituted 
into equation (E18) gives 
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(a) Stream surface over four stages of multistage 

tu rbomach 1 ne. 




(b) Intersection of 
stream surface with 
plane normal to axis. 


Ic) Intersection of stream 
surface with axial plane. 



Figure I. - Stream surface over four stages of multistage 
turbomachine and Its Intersection with planes normal to 
and containing axis of machine. 
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(a) Nontapered passage. 
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(b) Tapered passage. 



Figure 2. — Stations between blade rows 
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Figure 3. - Stations i and j short distance apart. 
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(b) Radial displacement across rotor. 


Figure 4. - Symmetrical-velocity-diagram and constant-total- 

enthalpy compressor. 
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(c) Variation of axial velocities. 


Figure 4. - Continued. Synraetricaf-velocity-diagran and constant-total -enthal py 

compressor. 
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r,/r t 


(ej Variation of air angles. 

Figure 4. - Continued. Symmet r i ca 1-vel oc i ty-d iag ram and constant- 

total-enthalpy compressor. 
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(c) Variation of axial velocities. 




Figure 5. - Continued. Free-vortex and constant-total-enthalpy compressor. 
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(e) Variation of Mach number relative to rotor blades. 

Figure 5. - Continued. Free-vortex and constant-total -enthal py 

compressor. 













CO 


> 

o 

> 





nthalpy compressor 


»D 


G/G 


NACA TN No 


95 



(a) Distribution of specific mass flow. 
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(b) Radial displacement across rotor. 



Figure 6. - Free-vortex and constant-total-enthalpy turbine 
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(e) Variation of'Mach number in front of rotor. 


Figure 6. - Continued. Free-vortex and constant-total-enthalpy 

tu rbi ne. 









